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1.  INTRODUCTION  AiiD  SUMMARY 

A major  mission  of  the  U.S.  Coast  Guard  is  the  task  of  providing  and 
maintaining  Maritime  Aids  to  Navigation.  These  aids  are  located  on  and 
near  the  coast  line  and  inland  waters  of  the  United  States  and  its  pos- 
sessions. A specific  type  of  navigational  aid  (flashing  lamp  cn  a buoy) 
offers  a possible  application  for  low-cost  solar  array/battery  power 
systems . 

Presently,  flashing  lamp  buoys  are  powered  by  zinc-air  batteries  which 
require  periodic  replacement  ( 1-2  year  life).  Development  of  a solar 
array/battery  power  system  with  unattended  operation  during  a six-year 
life,  shows  considerable  promise  for  reducing  buoy  maintenance  costs. 

Buoys  are  exposed  to  a wide  variety  of  environments.  The  range  of  these 
environments  are  typified  by  weather  conditions  in  locations  at  Fairbanks, 
Alaska,  and  Florida  coastal  waters.  Operational  and  design  problems 
are  compounded  when  it  is  realized  that  solar-powered  buoys  must  operate 
with  a wide  range  of  solar  insolation  in  addition  to  functioning  under  a 
variety  of  weather  conditions. 

The  Coast  Guard  has  a requirement  fcr  a computer  program  capable  of  per- 
forming all  the  calculations  necessary  to  understand  the  overall  charac- 
teristics of  solar  array/battery  power  systems  for  application  to  Aids  to 
Navigation.  This  report  contains  a description  of  the  Design  Synthesis/ 
Performance  Analysis  (DSPA)  Computer  Program  developed  by  the  Jet 
Propulsion  Laboratory  in  response  to  the  Coast  Guard  requirements. 

The  description  of  the  DSPA  Computer  Program  is  presented  as  indicated 
below: 

• Segment  A - Program  Summary. 

• Segment  B - Program  Documentation. 

Volume  I - Software  Requirements  Document. 

Volume  II  - Users'  Manual. 

Volume  III  - Programmer's  Manual. 
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.'he  DSPA  Computer  Program  will  enable  the  Coast  Guard  to: 

• Design  solar  array /battery  power  systems  for  aids  to  navigation 
having  any  specific  flashing  lamp  load  and  in  any  geographical 

, location. 

*•  . Simulate  the  performance  of  solar  array/battery  power  systems  in 

operation  on  navigational  aids  under  spiecific  flashing  lamp  loads 
and  environmental  conditions. 

• Provide  specifications  for  viable  alternative  designs  which  may  be 
used  to  integrate  solar  array/battery  power  systems  into  the  Coast 
Guard  inventory. 

The  specifications  generated  by  the  DSPA  computer  program  may  also  be  used 
by  the  Coast  Guard  as  a basis  for  a complete  economic  analysis  of  alterna- 
tive power  system  design. 

The  basis  for  selecting  the  computer  program  that  will  meet  the  Coast 
Guard  requirements  was  determined  by  an  examination  of  the  uses  of 
power  system  programs.  Generally,  power  system  computer  programs  are 
used  for  two  types  of  system  problems,  i.e.,: 

• Design  Synthesis. 

• Performance  Analysis. 

During  design  synthesis,  load  and  environmental  profiles  set  the  require- 
ments for  the  power  system.  Based  on  these  requirements  and  on  the 
electrical  characteristics  of  the  system  equipment,  the  computer  program 
determines  the  electrical  size  (volts,  amperes,  watts,  ampere-hours, 
watt-hours)  and  the  physical  characteristics  (weight,  area,  and  cost) 
of  the  power  system.  A design  synthesis  computer  program  can  also 
provide  an  easy  means  of  determining  sensitivities  of  the  power  system 
physical  characteristics  to  variations  in  the  many  equipment  parameters  - 
in  order  to  establish  the  rationale  for  "optimum"  power  system  designs. 

During  performance  analysis  calculations,  the  arrangements,  electrical 
size,  and  physical  characteristics  of  power  system  equipment  are  already 
known.  The  objective,  in  this  case,  is  to  determine  the  response  (or 
operational  characteristics)  of  the  power  system  equipment  to  a given 
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stimulus  (load  and  environmental  profiles).  The  performance  analysis 
computer  program  provides  a means  of  "measuring"  and  adjusting  steady 
state  performance  (and  transient  performance  if  the  selected  profile 
time  increments  are  small  enough)  of  proposed  or  existing  power  systems 
in  order  to  establish  designs  which  meet  system  requirements. 

The  determination  was  made  that  Coast  Guard  objectives  and  requirements  ! 

would  be  best  attained  if  the  computer  program  combined  the  features  of 
both  the  design  synthesis  and  performance  analysis  processes.  7 partic- 
ularly critical  requirement  for  determining  the  actual  daily  remaining 
battery  capacity  during  a one  year  period  was  particularly  amenable  to 
solution  by  the  performance  analysis  methodology. 

The  design  synthesis  methodology  allows  the  program  user  to  select  a lamp 

and  flasher  combination  as  well  as  a particular  power  system  arrangement  \ 

(out  of  a limited  set).  The  user  must  also  supply  information  on  daily 
cloud  cover  for  periods  of  one  week  up  to  fifty  two  (52)  weeks.  Based 
on  this  information,  the  computer  program  determines  the  hourly  solar 
insolation  as  well  as  the  day/night  load  durations  needed  for  a power 
load  profile.  (An  alternative  form  of  calculation  utilizes  solar  insola- 
tion data  taken  from  National  Oceanic  and  Atmospheric  (NOAA)  solar 
radiation  tapes).  The  power  load  profile,  after  modification  using 
battery  charge-efficiency  and  a number  of  power  system  cabling  and  diode 
losses,  is  used  in  the  load  profile  analysis.  The  object  of  a load 
profile  analysis  is  to  determine  the  electrical  size  of  a balanced* 
power  source  as  well  as  the  minimum  theoretical  electrical  size  of  a 
battery.  Once  this  information  is  obtained,  it  is  a fairly  simple  and 
straightforward  process  to  determine  the  electrical  size  of  the  remain- 
ing items  of  equipment  since  the  design  synthesis  methodology  is  a 
linear  step-by-step  procedure. 


^Balanced  power  source:  - power  source  which  provides  just  enough  energy  to 
the  batteries  (during  recharge  periods)  to  offset  or  balance  the  energy  loss 
sustained  by  the  batteries  during  periods  of  discharge. 
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The  methodology  used  in  performance  analysis*  is  to  obtain  the  raw  power 
bus  operating  point  for  each  selected  time  increment  during  a mission. 
Once  this  operating  point  is  obtained,  all  power  system  operational  char- 
acteristics are  easily  determinable. 

To  obtain  the  operating  point,  the  power  system  is  divided  into  groups 
about  the  operating  point  node.  These  groups  and  the  equipment  in  each 
group  are : 

Power  Source  Group: 

• Solar  Array. 

• Shunt  Limiter. 

• Solar  Array  Isolation  Diodes. 

• Solar  Array  Cable. 

Energy  Storage  Group: 

• Batteries. 

• Battery  Chargers. 

• Battery  Isolation  Diodes. 

• Battery  Cables. 

Power  Conditioning  and  Distribution  Group: 

• User  Loads. 

• Load  Cable. 

Current-voltage  (I-V)  characteristics  of  each  equipment  group  are  deter- 
mined. The  Power  Conditioning  and  Distribution  Group  characteristics  are 
then  deducted  from  those  of  the  Power  Source  Group  to  obtain  a Difference 
Curve.  The  intersection  of  the  Difference  Curve  with  the  Energy  Storage 
Group  characteristics  (as  determined  by  the  computer  program)  is  an 
estimate  of  the  voltage  and  current  on  the  raw  power  bus  (operating 
point).  During  the  time  increment  under  examination,  the  battery  is  in 
a state  of  charge,  discharge  or  open-circuit,  depending  on  the  value 
and  sign  of  the  operating  point  current.  After  the  operating  point  is 
obtained,  the  operational  characteristics  of  the  power  system  equipment 
are  determined.  The  computer  program  is  then  ready  for  examination  of 
power  system  response  for  the  next  time  increment.  The  process  is 
repeated  until  the  end  of  the  mission  of  interest. 


*This  portion  of  the  DSPA  Computer  Program  is  based  upon  a similar  program 
developed  for  the  analysis  of  planetary  spacecraft  power  systems  by  the  JPL 
and  for  the  National  Aeronautics  and  Space  Administration. 
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The  computed  results  of  both  design  synthesis  and  performance  analysis 
are  printed  out  in  easily  readable  and  comprehensive  tabular  formats. 
Additionally,  an  optional  graphical  summary  of  selected  performance  analy- 
sis output  data  is  available  as  an  aid  in  determining  performance  trends. 

A selection  of  key  data  used  to  represent  the  performance  characteristics 
of  the  equipment  was  compiled  and  also  presented.  This  program  data  is 
conveniently  placed  in  the  three  categories  representative  of  the  equip- 
ment functions,  i.e.,: 

• Power  Source  Group. 

• Energy  Storage  Group. 

• Power  Conditioning  and  Distribution  Group. 

The  final  phase  of  DSPA  Computer  Program  development  was  program  verifica- 
tion. This  phase  was  necessary  to  provide  a check  on  the  validity  of  the 
performance  analysis  methodology  used  in  the  program  by  comparing  the 
predicted  operational  characteristics  of  a known  power  system  with  the 
actual  operational  performance  of  that  same  power  system.  During  this 
phase,  data  on  the  critical  parameters  of  five  power  systems,  undergoing 
test  at  the  USCG  Research  and  Development  Center,  Groton,  Conn.,  were 
supplied  by  the  Coast  Guard. 

The  DSPA  Computer  Program  accurately  predicted  the  change  in  battery  capa- 
city over  a one  year  period  for  the  two  power  systems  considered  success- 
fully tested  by  the  Coast  Guard.  The  other  three  power  systems  exper- 
ienced erratic  solar  array  performance  and  were  considered  failures  by 
the  Coast  Guard.  In  spite  of  these  problems,  the  DSPA  fairly  well 
predicted  the  battery  capacity  over  a one  year  period  for  two  of  the  three 
"failed"  power  systems. 
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2.  COMPUTER  PROGRAM  DESCRIPTION 


The  features  of  the  DSPA  Computer  Program  are  discussed  in  this  section 
of  the  program  summary.  The  key  elements  of  this  discussion  include 
specific  program  requirements,  computer  program  rationale,  DSPA  Computer 
Program,  design  synthesis  and  performance  analysis.  These  elements  are 
presented  below. 


2.1  Specific  Program  Requirements 


The  major  requirements  of  the  DSPA  Computer  program,  as  discussed  in  the 

Introduction  were  to  provide  a capability  to: 

a.  Design  solar  array/battery  power  systems  for  aids  to  navigation  having 
any  specific  flashing  lamp  load  and  at  any  geographical  location. 

b.  Simulate  the  performance  of  solar  array /battery  power  systems  in 
operation  on  Navigational  Aids  using  specific  flashing  lamp  loads, 
at  any  geographical  location  and  under  specific  environmental 
conditions . 

c.  Determine  specification  for  viable  alternative  designs  which  may  be 
used  to  integrate  solar  array/battery  power  systems  into  the  Coast 
Guard  inventory. 

Additionally,  there  were  specific  requirements  that  the  DSPA  Computer 

Program  have  the  following  features: 

d.  Program  flexibility  to  accept  data  on  solar  insolation,  wind  and 
wave  motion.  It  was  intended  that  the  program  be  written  with 
enough  flexibility  so  that  power  system  design  analysis  can  be 
determined  for  systems  using  solar  insolation,  wind  or  wave  motion 
as  the  primary  energy  sources.  The  initial  version  of  the  computer 
program,  however,  was  only  required  to  utilize  solar  array  power 
sources  with  the  capability  of  easily  adding  subprograms  which 
utilize  the  characteristics  of  wind-powered  or  wave  motion  power 
generators . 

e.  Ability  to  generate  solar  insolation  values  at  any  latitude  and 
longitude. 

f.  Capability  of  accepting  input  parameters  such  as  temperature,  hours 
of  sunlight,  load  characteristics,  and  cost  factors.  JPL  was 
further  required  to  identify  any  additional  input  parameters  needed 
for  program  operation.  As  an  example,  the  use  of  this  computer 
program  for  both  buoys  and  shore  stations  was  desired.  In  the  case 
of  buoys,  the  effects  of  wind/wave  motion  on  array  power  output  (and 
other  characteristics)  were  to  be  ignored.  The  assumption  was  to 
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be  made  that  all  buoy  solar  arrays  are  horizontal.  Shore  station 
arrays,  on  the  other  hand,  would  be  southfacing  and  tilted  at  some 
"optimum"  angle  (e.g.,  tilt  angle  't  latitude),  lienee,  the  computer 
program  was  coded  to  accept  any  preselected  solar  array  tilt  angle. 

Program  modular  design,  wherein  the  characteristics  f each  eq;ip- 
ment  module  could  be  readily  included,  deleted,  r changed  by 

the  computer  program  operator.  Equipment  modules  include  solar 
cell  panels,  batteries,  battery  chargers,  shunts,  etc.  The  program 
operator,  for  example,  could  select  the  type  of  shunt  desired, 
i.e.,  none,  ordinary  zener  diode,  temperature-compensated  zener 
diode  or  active  type.  Characteristics  of  all  of  these  shunt  types 
are  part  of  the  normal  program  data  which  may  be  modified  by  the 
program  operator  as  new  characteristics  are  determined. 

Automatic  design  operations  with  the  capability  of  allowing  the 
program  operator  to  specify  variable  operating  parameters  and 
design  criteria.  Examples  of  design  criteria  over  which  the 
operator  must  have  control  include: 

1)  Minimum  remaining  acceptable  battery  capacity. 

2)  Inclusion  of  restrictions  on  solar  array  and  battery  sizing. 

3)  Inclusion  of  safety  factors. 

1+)  Selection  of  system  elements. 

Perform  calculations  which  take  into  consideration  all  parameters 
affecting  useable  stored  battery  capacity.  Examples  of  these 
parameters,  listed  below,  are  typically  obtained  from  analyses  of 
power  system  operational  characteristics. 

1)  Solar  Array  current-voltage  characteristics  as  a function  of 
solar  insolation  magnitude  and  cell  temperature. 

2)  Voltage  Regulator:  Voltage  cut-off  characteristic  as  a func- 
tion of  input  voltage  and  temperature. 

3)  Storage  Battery:  Lifetime  degradation  characteristics  (cycle 
life)  of  the  battery  as  a function  of  depth  of  discharge. 

Capability  of  storing  or  generating  the  total  range  of  operating 
characteristics  of  all  possible  system  elenents. 

Ability  to  accept  various  restrictions  on  power  system  operations. 
Examples  of  these  restrictions  are: 

l)  Minimum  battery  temperature  as  a function  of  remaining  capacity. 
Information  on  freezing  temperatures  are  generally  obtained 
from  data  of  the  type  shown  in  Figure  2-1. 
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(a)  PERCENT  CAPACITY  REMAINING 


(b)  FREEZING  TEMPERATURE 


FIGURE  2-1.  BATTERY  ELECTROLYTE  FREEZING  TEMPERATURE  AS  FUNCTION 
OF  REMAINING  CAPACITY 

2)  Minimum  allowable  remaining  battery  capacity. 

3)  Maximum  time  allowable  with  the  battery  capacity  below  a given 
value . 

l.  Program  Outputs.  The  primary  program  output  must  be  a table  which 
shows  actual  remaining  battery  capacity  for  each  day  over  a one  to 
two  year  period.  This  is  the  central  and  critical  requirement  for 
the  entire  computer  program.  The  computer  program  must  be  capable 
of  dealing  with  a one  year  load/environmental/solar  insolation 
profile  for  a reasonable  running  cost. 

m.  Program  Language.  The  computer  program  must  be  written  in 
FORTRAN-V  for  use  on  a UNIVAC-1108  Computer.  Such  a computer  is 
available  to  personnel  from  the  U.S.  Coast  Guard  R&D  Center,  Avery 
Point,  Conn.,  and  is  located  at  the  Naval  Underwater  Systems 
Center  (NUSC),  New  London,  Conn. 


2.2  Computer  Program  Rationale 


The  basis  for  selecting  the  proper  computer  program  that  will  meet  the 
Coast  Guard  requirements  was  determined  by  an  examination  of  the  uses 
of  power  system  programs.  Generally,  power  system  computer  programs  are 
used  for  two  types  of  system  problems,  i.e.,: 

• Design  Synthesis. 

• Performance  Analysis 
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In  design  synthesis,  load  and  environmental  profiles  set  the  requirements 
for  the  power  system.  Based  on  these  requirements  and  on  the  electrical 
characteristics  of  the  system  equipment,  the  computer  program  determines 
the  electrical  size  (volts,  amperes,  watts,  ampere-hours,  and  watt-hours) 
and  the  physical  characteristics  (weight,  area,  cost)  of  the  power  system. 

A design  synthesis  computer  program  can  also  provide  an  easy  means  of 
determining  sensitivities  of  the  power  system  physical  characteristics  to 
variations  in  the  many  of  the  equipment  parameters-in  order  to  establish 
the  rationale  for  "optimum"  power  system  designs.  Some  of  the  more  flex- 
ible design  synthesis  programs  may  even  select  the  "optimum"  arrangments 
of  power  system  equipment,  rather  than  leaving  the  choice  of  such 
arrangements  to  the  program  operator. 

In  performance  analysis  calculations,  the  arrangements,  electrical  size 
and  physical  characteristics  of  power  system  equipment  are  already  known. 
The  objective,  in  this  case,  is  to  determine  the  response  (or  operational 
characteristics)  of  the  power  system  equipment  to  a given  stimulus  (load 
and  environmental  profiles).  The  performance  analysis  computer  program 
provides  a means  of  "measuring"  and  adjusting  steady-state  performance 
(and  transient  performance  if  the  selected  profile  time  increments  are 
small  enough)  of  proposed  or  existing  power  systems  in  order  to  establish 
designs  which  meet  system  requirements. 

A number  of  methods  for  performance  analysis  of  circuits  and  modules  have 
been  developed  in  the  past.  Usually,  these  take  the  form  of  nodal  analyses 
in  which  the  transfer  function  of  a circuit,  part,  or  component  is  analy- 
zed in  a complex  network  whose  branch  interconnection  points  or  nodes  are 
reference  points.  A voltage  or  current  is  assumed  for  each  node,  the 
entire  circuit  solved  by  solution  of  simultaneous  equations,  and  the 
values  of  voltage  and  current  at  each  node  compared  with  the  assumed 
values.  If  they  are  not  the  same,  the  computed  values  replace  the  initial 
assumptions,  and  the  solution  is  repeated.  After  a number  of  iterations, 
the  system  eventually  converges,  and  the  final  solution  is  reached.  If 
the  system  consists  of  more  than  a very  few  nodes,  the  problem  is  solvable 
only  in  a computer  and  with  large  numbers  of  nodes.  Even  the  computer 
programs  can  become  too  costly  to  use  because  of  the  increasing  number  of 
iterations  required  to  reach  a final  convergence.  When  these  problems  are 


2-U 


501+0-27 


added  to  the  necessity  for  possible  thermal  solutions  at  each  moment  in 
time,  use  of  nodal  analysis  programs  for  describing  the  performance  of 
power  subsystems  containing  batteries  presents  a discouraging  picture. 

The  approach  herein  to  Performance  Analysis,  discussed  in  Section  2.5 
below,  represents  an  achievable  compromise  between  the  desirable  complete 
electrical  analysis  given  by  the  nodal  analysis  type  of  program  and  the 
excessive  cost  and  complexity  of  detailed  nodal  analysis.  The  compromise 
approach  (to  Performance  Analysis)  uses  a graphical  analysis  method  which 
is  implemented  for  the  computer  and  which  also  lends  itself  to  more 
accurate  modeling  of  battery  current-voltage  characteristics. 

The  determination  was  made  that  Coast  Guard  objectives  and  requirements 
would  be  best  attained  if  the  computer  program  combined  the  features  of 
both  the  design  synthesis  and  the  compromise  approach  to  performance 
analysis  processes.  The  particularly  critical  requirement  for  determining 
the  actual  daily  remaining  battery  capability  during  a one  year  period  was 
particularly  amendable  to  solution  by  the  "compromise"  performance 
analysis  methodology. 

Accordingly,  the  computer  program  actually  developed  was  one  which  com- 
bines the  elements  of  design  synthesis  and  performance  analysis  and  was 
named  the  Design  Synthesis  Performance  Analysis  (DSPA)  Computer  Program. 

2.3  DSPA  Computer  Program 

The  generalized  arrangement  of  the  solar  array/battery  power  system  used 
as  the  basis  of  the  DSPC  Computer  Program  is  shown  in  Figure  2-2.  Basic- 
ally, the  system  contains  the  following  major  elements: 

• Solar  array. 

• Shunt  Limiter. 

• Battery  charger. 

• Battery. 

• User  Loads. 
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FIGURE  2-2.  POWER  SYSTEM  BLOCK  DIAGRAM 


Solar  arrays  are  the  primary  power  source,  providing  for  the  photovoltaic 
energy  conversion  of  sunlight  into  electrical  power.  The  array  must  be 
large  enough  to  provide  sufficient  power  for  the  user  loads,  power  system 
"housekeeping"  losses,  and  battery  charging  during  periods  of  solar 
insolation. 


A shunt  limiter  may  be  provided  to  limit  the  maximum  voltage  on  the  raw 
power  bus.  The  shunt  limiter  is  usually  designed  to  accept  the  full 
output  power  of  the  solar  array. 


Batteries  are  needed  to  provide  power  during  periods  when  the  user  loads 
require  more  power  than  is  available  from  the  solar  array  as  well  as  for 
providing  power  during  periods  of  solar  occultation.  The  minimum  voltage 
on  the  raw  power  bus  is  determined  by  the  minimum  battery  voltage. 


Battery  chargers  may  also  be  used  in  the  process  of  restoring  energy  to 
a discharged  battery.  The  function  of  the  charger  is  to  maintain  battery 
voltage  and  current  within  tolerable  levels  during  the  charging  process. 
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The  user  loads  consist  of  a flashing  lamp  (utilized  on  a buoy)  as  well 
as  a regulator  to  control  the  current  voltage  and  phasing  of  that  lamp. 

An  important  concept  in  power  system  operation  and  analysis  is  that  of  the 
the  operating  point,  i.e.,  the  nodal  point  in  the  power  circuit  about 
which  the  major  elements  of  the  power  system  are  grouped.  As  an  aid  in 
design  and  analysis,  the  system  elements  are  divided  into  three  major 
9 groups  about  the  operating  point  node. 

• Power  Source  Group: 

Solar  Array. 

Shunt  Limiter. 

Solar  Array  Isolation  Diodes. 

Solar  Array  Cable. 

• Energy  Storage  Group: 

Batteries . 

Battery  Chargers. 

Battery  Isolation  Diodes. 

Battery  Cables. 

• Power  Conditioning  and  Distribution  Group: 

User  Loads . 

Load  Cable. 


The  Energy  Storage  Group  consists  of  parallel  Energy  Storage  Units.  A 
single  Energy  Storage  Unit,  as  shown  in  Figure  2-2,  contains  a battery, 
battery  charger,  battery  isolation  diode  and  battery  cables. 


A functional  block  diagram  of  the  Design  Synthesis  Performance  Analysis 
(DSPA)  Computer  Program  is  shown  in  Figure  2-3.  The  DSPA  Computer  Pro- 
gram, at  the  primary  level,  utilizes  the  following  methodology: 

a.  The  program  operator  selects  the  power  system  arrangement  desired. 

b.  If  a design  synthesis  i_s  not  required,  then  the  program  operator 
must  provide  information  on  the  electrical  size  of  the  equipment 
selected. 

c.  If  a design  synthesis  is_  requested,  the  program  operator  must 
supply  information  on  the  parameters  used  in  determining  the  various 
profiles.  The  computer  program  then  calculates  the  load  and  environ- 
mental profiles  needed  in  their  entirety  for  a load  profile 
analysis.  Based  on  a profile  energy  balance  determined  as  part  of 
the  load  profile  analysis,  the  computer  estimates  the  electrical 
size  of  the  equipment  selected. 


L 
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FIGURE  2-3.  DESIGN  SYNTHESIS  PERFORMANCE  ANALYSIS  COMPUTER  PROGRAM 


d.  The  physical  characteristics  of  the  equipment  selected  is  then 
determined  by  the  computer,  based  on  the  equipment  electrical  size. 

e.  If  performance  analysis  is_  not  required,  the  data  calculated  by  the 
program,  along  with  the  input  data,  is  printed  in  the  appropriate 
output  data  format. 

f.  If  performance  analysis  requested,  the  program  operator  must 
again  supply  information  on  the  parameters  used  in  determining  the 
various  profiles.  The  computer  program  then  calculates  the  values 
of  the  load  and  environment  at  the  start  of  a selected  mission 
period.  These  stimuli  (load,  environment)  are  used  to  calculate 
the  response  (operational  characteristics)  of  the  equipment  at  that 
point  in  time.  The  process  is  repeated  for  selected  time  increments 
until  the  power  system  operational  characteristics  for  the  entire 
mission  period  have  been  determined.  This  information  is  then 
printed  in  the  appropriate  output  data  format. 


An  overview  of  the  input/output  elements  of  the  DSPA  Computer  Program 
is  shown  in  Figure  2-1*.  Input  data  requirements  fall  into  two  major 
groups,  i.e.,  Program  Parameter  Data  and  Environmental  Profiles. 


Program  Parameter  Data  contains  information  on  mission  requirements 
(latitude,  longitude,  duration,  etc)  and  equipment  characteristics. 
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FIGURE  2-4.  DSPA  COMPUTER  PROGRAM  INPUT/OUTPUT 

These  data  are  required  in  a "Namelist"  format  - for  ease  in  changing 
values  of  the  parameters  for  repetitive  computer  runs. 


Environmental  Profiles  are  obtained  from  two  sources.  The  program  user 
may  elect  to  prepare  his  own  cloud  cover  profile  (in  Free-Format)  for 
use  in  estimating  solar  insolation,  or  he  may  elect  to  use  the  actual 
solar  insolation  values  for  a given  location  directly  from  a "MERGE" 
tape  (described  below). 


After  DSPA  Computer  Program  calculations  are  completed,  the  output  will 
be  printed  in  a tabular  format.  The  program  user  may  also  elect  to 
receive  a graphical  output  of  the  tabulated  results. 


Solar  Insolation  and  other  environmental  information  is  available  on 
magnetic  tapes  from  the  National  Climatic  Center,  National  Oceanic  and 
Atmospheric  Administration  (NOAA)  at  Asheville,  North  Carolina.  The 
magnetic  tapes  of  particular  interest  are  those  for  Surface  Observations 
(TDF-14;  Reference  l)  and  Solar  Radiation  (DECK-280;  Reference  2). 
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The  information  recorded  on  TDF-lb  magnetic  tapes  is  given  in  Table  2-1. 
Data  recorded  on  Deck  280  magnetic  tapes  are  shown  in  Table  2-2. 


Since  these  tapes  contain  more  information  than  was  needed  as  input  for 
the  DSPA  Computer  Program,  an  independent  "MERGE"  program  was  developed 
to  provide  a simplified  tape  input.  The  "MERGE"  program  performs  the 
following  functions. 

• Reads  in  data  from  the  TDF-lU  and  Deck-280  tapes. 

• Combines  data  into  single,  sequential  random  access  files. 

• Annotates  empty  data  locations. 

• Prepares  a "MERGE"  tape  containing  temperature,  solar  insolation 
and  wind  velocity  as  a function  of  time  for  various  geographical 
locations . 


"MERGE"  tapes  have  been  prepared  for  the  ten  principal  coastal  locations 
(shown  in  Table  2-3)  and  five  supplemental  locations  (shown  in  Table  2-1+) 
selected  by  the  Coast  Guard. 


TABLE  2-1. -SURFACE  OBSERVATIONS,  TDF-ll*  INFORMATION  RECORDED 


Station  Humber 

Drizzle 

Blowing  Spray 

Year,  Month,  Day,  Hour 

Freezing  Drizzle 

Sea  Level 

Ceiling  Height 

Snow 

Pressure 

Sky  Condition: 

Snow  Pellets 

Dew  Point 

Clear 

Ice  Crystals 

Temperature 

Scattered 

Snow  Showers 

Wind  Direction 

Broken 

Snow  Grains 

Wind  Speed 

Overcast 

"Sleet 

Station  Pressure 

Partial  Obscuration 

Hail 

Dry  Bulb 

Obscuration 

"Small  Hail 

Temperature 

Visibility 

Fog 

Wet  Bulb 

Weather  and/or 

Ice  Fog 

Temperature 

Obstruction  to  Vision 

Ground  Fog 

Relative  Humidity 

Thunde rstorm 

Blowing  Dust 

Total  Sky  Cover 

Tornado 

Blowing  Sand 

Amount,  Type  and 

Squall 

Smoke  and/or 

Height  of  Cloud 

Rain 

Haze 

Layers 

Rain  Showers 

Dust 

Opaque  Sky 

Freezing  Rain 

Blowing  Snow 

Cover 

"Reported  as  Ice  Pellets, 

effective  1 Apr  70. 
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iABLJ ; 2-2. -SOLAR  RADIATION,  DECK-280  INFORMATION  RECORDED 


Station  Number 

Year,  Month,  Day,  Hour 

Solar  Radiation  — Hemispheric  (sum  of  direct  and  diffuse) 

Solar  Elevation 

Extra-Terrestrial  Radiation 

Sunshine 

Snow  Cover 

Solar  Week  (discontinued  1 Jan  6 3) 

Opaque  Sky  Cover  (discontinued  1 Jan  65) 

Solar  Hour 

Percent  of  Possible  Radiation  (discontinued  1 Jan  65) 

Visibility  (discontinued  1 Jan  6 5) 

WeaUiUT  hihI'AjY  Ht3STTU<il.l6ns  to  vision'Caiscontinued  1 Jan  65) 

Dry  Bulb  Temperature  °F  (discontinued  1 Jan  65) 

Dew  Point  Temperature  °F  (discontinued  1 Jan  65) 

Amount,  Type  and  Height  of  Cloud  Layers  (discontinued  1 Jan  65) 

2.1+  Design  Synthesis 


A detailed  picture  of  the  design  synthesis  methodology  is  shown  in  Fig- 
ure 2-5.  This  methodology  allows  the  programmer  user  to  select  a lamp  and 
flasher  combination  as  well  as  a particular  power  system  arrangement  (out 
of  a limited  set).  The  user  must  also  supply  information  on  daily  cloud 
cover  for  periods  of  one  week  up  to  fifty  two  (52)  weeks  or,  if  so  desired, 
hourly  solar  insolation  and  temperature  values  from  a "MERGE"  tape  for  the 
same  periods  of  time.  Based  on  this  information,  the  DSPA  Computer  Pro- 
gram determines  the  daily  periods  of  solar  occultation,  power  system  share- 
mode operation,  battery  charging  and  the  number  of  battery  operating  mode 
reversals  (needed  to  estimate  battery  charge/discharge  cycle  requirements). 
Using  the  information,  a power  load  profile  is  obtained.  The  power  load 
profile,  after  modification  using  battery  charge-efficiency  and  a number 
of  power  system  cabling  and  diode  losses,  is  used  in  the  load  profile 
analysis.  The  object  of  a load  profile  analysis  is  to  determine  the  elec- 
trical size  of  a balanced*  power  source  as  well  as  the  minimum  theoretical 
size  of  a battery.  Mathematical  formulations  for  estimating  the  size  of  a 


^Balanced  power  source:  -a  power  source  which  provides  just  enough  energy  to 
the  batteries  (during  recharge  periods)  to  offset  or_  balance  the  energy  loss 
sustained  by  the  batteries  during  periods  of  discharge. 
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TABLE  2-3 • -liOAA  SOLAR / WEATHEK  TAPES 


^ .Pi’ a. nc ip cil  ijCCu,uxw.i j / 


Station 

Number 

Station 

Name 

NOAA 

Tape 

Data  Periods  Covered 

1U607 

Caribou,  Maine 

DK-280 

Jan  1955  - Dec  196U 

11+607 

Caribou,  Maine 

tdf-il 

Jan  1955  - Dec  1961* 

ll*753 

Blue  Hill/Milton,  Mass. 

DK-280 

Jan  1955  - Dec  1961* 

11*739 

Boston,  Mass. 

TDF-ll* 

Jan  1955  - Dec  1961* 

91*706 

New  York,  N.Y.,  (Central  Park) 

DK-280 | 

Jan  1953  - Dec  1958 

91*706 

New  York,  N.Y.,  (Central  Park) 

DK-280 1 

Jan  1963  - Dec  1966 

11*732 

La  Guardia  Field,  N.Y. 

TDF-lL  | 

Jan  1953  - Dec  1958 

11*732 

La  Guardia  Field,  N.Y. 

TDF-lb  1 

Jan  1961  - Dec  1961+ 

11*81*7 

Sault  Ste.  Marie,  Mich. 

DK-280 

July  1952  - Aug  1958 

11*81*7 

Sault  Ste.  Marie,  Mich. 

TDF-ll* 

Jan  1952  - Dec  1961 

1371*5 

Hatteras , N.C. 

DK-280 ) 

Jan  1955  - Mar  1957 

93729 

Cape  Hatteras,  N.C. 

DK-280  l 

Mar  1957  - Dec  1961+ 

1371*5 

Hatteras,  N.C. 

TDF-ll*  1 

Jan  1955  - Mar  1957 

93729 

Cape  Hatteras,  N.C. 

TDF-ll*  I 

Mar  1957  - Dec  I96U 

12839 

Miami,  Florida 

DK-280 

Jan  1955  - Dec  1961+ 

12839 

Miami,  Florida 

TDF-lL 

Jan  1955  - Dec  19 61* 

12919 

Brownsville,  Texas 

DK-280 | 

Jan  1953  - Dec  1955 

12919 

Brownsville,  Texas 

DK-280  1 

Jan  1959  - Dec  1965 

12919 

Brownsville,  Texas 

TDF-lL  | 

Jan  1953  - Dec  1956 

12919 

Brownsville,  Texas 

TDF-ll*  1 

Jan  1959  - Dec  196I+ 

23171* 

Los  Angeles,  Calif. 

DK-280 

Jan  1962  - Dec  1966 

23171* 

Los  Angeles,  Calif. 

TDF-ll* 

Jan  1955  - Dec  196I+ 

21*233 

Seattle,  Wash. 

DK-280 

Jan  1955  - Dec  I96U 

2U233 

Seattle,  Wash. 

TDF-ll* 

Jan  1955  - Dec  196I+ 

26615 

Bethel,  Alaska 

DK-280  I 

July  1952  - Oct  1952 

26615 

Bethel,  Alaska 

DK-280 1 

Dec  1956  - Apr  1957 

26615 

Bethel,  Alaska 

TDF-ll* 

Jan  1952  - Dec  1961 

balanced  power  source  are  given  in  Appendix  A.  The  equations  in  Appendix  A 
form  the  basis  for  the  key  algorithms  used  in  design  synthesis. 


Once  the  electrical  size  of  a balanced  power  source  is  obtained,  it  is  a 
fairly  simple  and  straight  forward  process  to  determine  the  electrical 
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TABLE  2-U.-N0AA  SOLAR/ WEATHER  TAPES 


(Supplemental  Locations) 


Station 

Number 

Station 

Name 

NOAA 

Tape 

Data  Periods  Covered 

13983 

Columbia,  MO. 

DK-280 

Jan  1955  - Dec  1964 

13983 

Columbia,  MO 

TDF-14 

Jan  1955  - Dec  1964 

23154 

Ely,  NV. 

DK-280 

Jan  1955  - Dec  1964 

23154 

Ely,  NV. 

TDF-14 

Jan  1955  - Dec  1964 

24225 

Medford,  OR. 

DK-280 

Jan  1955  - Dec  1964 

24225 

Medford,  OR. 

TDF-14 

Jan  1955  - Dec  1964 

93193 

Fresno,  CA. 

DK-280 

Jan  1955  - Dec  1964 

93193 

Fresno,  CA 

TDF-14 

Jan  1955  - Dec  1964 

93722 

Silver  Hill,  MD. 

DK-280 | 

Jan  1955  - Dec  i960 

93734 

Sterling,  VA. 

DK-280 ) 

Dec  i960  - Dec  1964 

13743 

Wash.  D.C. 

(Natl.  Airport) 

TDF-14 

Jan  1955  - Dec  1964 

size  of  the  remaining  items  of  equipment  since  the  design  synthesis 
methodology  is  a linear  step-by-step  procedure. 


The  major  elements  (subroutines)  of  the  design  synthesis  subprogram  are 
shown  in  Figure  2-6.  Functions  performed  by  each  of  the  major  subroutines 
are  also  shown  in  this  figure. 


Examples  of  tabular  output  from  design  synthesis  calculation  are  shown  in 
Figures  2-7  to  2-10.  The  output  includes  separate  tabular  parameters  for: 

• Design  specifications  and  costs. 

• Engineering  Design  Data. 

• Battery  Performance  Data. 

• Power  Load  Profile  Analysis. 


The  first  two  of  the  presentation  shows  the  costs,  specifications  and 
major  physical  characteristics  of  the  designed  power  system.  The  last 
two  presentations  cover  key  information  obtained  during  tne  design  process. 


2-13 


FIGURE  2-5.  POWER  SYSTEM  DESIGN  SYNTHESIS 
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MAIN  DESIGN  SYNTHESIS  DRIVER  PROGRAM 


DS-BATEFC  COMPUTES  BATTERY  EFFICIENCY 


DS-CDSI  COMPUTES  CLEAR  DAY  SOLAR  INSOLATION 


COMPUTES  BATTERY  CURRENT  AND  VOLTAGE  CHARACTERISTICS 


COMPUTES  POWER  CONDITIONING  AND  DISTRIBUTION 
GROUP  CURRENT  AND  VOLTAGE  CHARACTERISTICS 


COMPUTES  SOLAR  ARRAY  CURRENT  AND  VOLTAGE  CHARACTERISTICS 


COMPUTES  TERMINATOR  CHARACTERISTICS 
(SUNRISE  AND  SUNSET  TIMES) 


FIGURE  2-6.  DESIGN  SYNTHESIS  SUBPROGRAM  ELEMENTS 
2.5  Performance  Analysis 


A detailed  picture  of  the  performance  analysis  methodology  is  shown  in 
Figure  2-11.  This  portion  of  the  DSPA  Computer  Program  was  based  upon  a 
similar  program  developed  by  JPL  for  the  analysis  of  planetary  spacecraft 
power  systems  (Reference  3). 


The  analysis  of  terrestrial  power  system  performance  is  a complex  task. 

This  complexity  stems,  in  part,  from  the  interaction  of  the  power  system 
with  the  ambient  environment.  Batteries,  for  example,  have  electrical 
characteristics  which  show  significant  variation  with  temperature.  The 
same  is  true  of  solar  arrays.  Hence,  at  a minimum,  any  reasonably  accurate 
simulation  of  power  system  operation  must  include  ambient  and  equipment 
temperatures  as  an  independent  parameter.  An  overview  of  the  complexity 
of  the  procedures  necessary  to  analyze  power  system  performance  is  indi- 
cated in  Figure  2-11.  This  figure  shows  the  flow  of  lata  during  the 
calculation  as  well  as  the  iterative  nature  of  the  procedures. 
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FIGURE  2-8.  ENGINEERING  DESIGN  DATA 
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FIGURE  2-9.  BATTERY  PERFORMANCE  DATA 
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FIGURE  2-1C . POWER  LOAD  PROFILE  ANALYSIS  (Sheet  1 of  2) 
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FIGURE  2-10.  POWER  LOAD  PROFILE  ANALYSIS  (Sheet  2 of  2) 
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The  performance  analysis  methodology  uses  a graphical  analysis  method, 
implemented  for  the  computer  to  obtain  the  power  system  operating  point 
for  each  selected  time  increment  during  a mission. 


Once  this  operating  point  is  obtained,  all  power  system  operational  charac- 
teristics are  easily  determinable.  To  obtain  the  operating  point,  the 
power  system  is  divided  into  groups  about  the  operating  point  (supra). 

These  groups  are: 

• Power  Source  Group  (PSG). 

• Energy  Storage  Group  (ESG). 

• Power  Conditioning  and  Distribution  Group  (PCDG). 

Current-voltage  (I-V)  characteristics  of  each  equipment  group  are  deter- 
mined. The  Power  Concitioning  and  Distribution  Group  characteristics  are 
then  deducted  from  those  of  the  Power  Source  Group  to  obtain  a Difference 
Curve.  The  intersection  of  the  Difference  Curve  with  the  Energy  Storage 
Group  characteristics  (as  determined  by  the  computer  program)  is  an 
estimate  of  the  voltage  and  current  on  the  raw  power  bus  (operating  point). 
During  the  time  increment  under  examination,  the  battery  is  in  a state 
of  charge,  discharge  or  open-circuit,  depending  on  the  value  and  sign  of 
the  operating  point  current.  After  the  operating  point  is  obtained,  the 
operational  characteristics  of  the  power  system  equipment  are  determined. 
The  computer  program  is  then  ready  for  examination  of  power  system 
response  for  the  next  time  increment.  The  process  is  repeated  until  the 
end  of  the  mission  of  interest. 

Power  Source  Group  characteristics  are  determined  as  shown  in  Figure  2-12. 
As  a starting  point,  solar  array  characteristics  are  obtained.  These 
characteristics  are  based  on  the  analytical  procedures  developed  for  the 
JPL  Planetary  Spacecraft  Power  Systems  Program  modified  and  incorporated 
into  the  DSPA  computer  program. 

A major  parameter  needed  in  determining  solar  array  characteristics  is 
the  level  of  solar  insolation.  Solar  insolation  nay  be  obtained  directly 
from  a "l-IERGE"  tape  (supra)  or  estimated  by  the  DSPA  computer  program. 
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The  methodology  used  in  making  solar  insolation  estimates  is  discussed  in 
Appendix  B.  The  techniques  were  based  upon  formulations  and  algorithms 
developed  by  the  American  Society  of  Heating,  Refrigerating  and  Air  Condi- 
tioning Engineers  (Reference  4).  Solar  array  characteristics  are  further 
determined  by  single  solar  cell  characteristics  and  degradation  as  well  as 
the  series  - parallel  arrangement  of  cells  in  each  module. 

Shunt  Limiter  characteristics  are  determined  by  the  type  of  Limiter 
selected,  the  desired  raw  power  bus  voltage  cutoff  and  the  limiter 
temperature . 

The  modified  solar  array  characteristic,  as  shown  in  Figure  2-12a  is 
obtained  by  subtracting  the  shunt  limiter  current  from  the  solar  array 
current  at  each  point  between  zero  volts  and  array  open  circuit  voltage. 

Power  Source  Group  characteristics,  shown  in  Figure  2-12b,  are  obtained  by 
subtracting  the  isolation  diode  and  cable  voltage  drops  from  the  modified 
solar  array  voltage  at  each  point  between  zero  amperes  and  modified  solar 
array  short  circuit  current. 

The  methodology  of  determining  Power  Conditioning  and  Distribution  Group 
characteristics  is  discussed  in  Appendix  C.  Using  these  characteristics, 
the  interaction  between  the  Power  Source  Group  and  the  Power  Conditioning 
and  Distribution  Group  is  shown  in  Figure  2-13. 

A load  line  difference  curve  is  used  to  determine  the  interaction  of  the 

i 

Power  Source,  Energy  Storage  and  PC&D  Groups.  The  difference  curve,  as 
indicated  in  Figure  2-13b,  is  obtained  by  subtracting  the  PC&D  Group  cur- 
rent from  the  Power  Source  Group  current  at  each  voltage  level.  By  oper- 
ating along  positive  values  of  the  difference  curve,  the  solar  array  has 
excess  power  which  can  be  used  to  charge  the  batteries.  Operating  on  neg- 
ative values  of  the  difference  curve  will  result  in  battery  discharge. 

Techniques  for  estimating  the  current-voltage  characteristics  of  the 
Energy  Storage  Group  are  giver;  in  Appendix  D.  Using  these  characteristics, 
the  interaction  between  the  Energy  Storage  Group  and  the  Load  Difference 
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Curve  during  Load  Line  Analysis  is  shown  in  Figure  2-1*4.  Load  Line  Analy- 
sis is  performed  by  determining  the  point  of  intersection  (power  system 
operating  point)  of  the  Energy  Storage  Group  Characteristic  with  the  Load 
Difference  Curve. 

The  interactions  in  Figure  2-1*4  also  indicate  the  change  in  power  system 
operation  with  battery  state  of  charge.  These  several  operating  points 
include : 

• Battery  Charging  - (Operating  Points  ID,  2D). 

• Battery  Floating  - (Operating  Point  3D). 

• Share  Mode  (Battery  Discharging  and  sharing  load  with  Solar  Array) 

- (Operating  Point  *4D). 

• Battery  Discharge  - (Operating  Points  5D,  6D,  IN-6ll). 

The  major  elements  (subroutines)  of  the  performance  analysis  subprogram 
are  shown  in  Figure  2-15.  Functions  performed  by  each  of  the  major  sub- 
routines are  also  shown  in  this  figure. 

Examples  of  tabular  output  from  performance  analysis  calculations  are 
shown  in  Figures  2-l6  to  2-18.  This  output  includes  separate  tabular 
presentations  for: 

• Unregulated  Bus  Summary  Data. 

• Power  Source  Group  Summary  Data. 

• Energy  Storage  Unit  Summary  Data. 

These  tabulation  show  the  predicted  performance  of  the  major  power  system 
elements  for  selected  time  increments  over  a selected  period  of  the 
mission. 

One  form  of  graphical  output  is  shown  in  Figures  2-19  to  2-21.  These 
forms  of  output  illustrate  the  factors  involved  in  load  line  analysis  for 
selected  times  during  the  mission.  The  particular  times  shown  in  these 
figures  are  at  *4:00  a.m. , 8:00  a.m.  and  12:00  noon  on  January  1,  1975. 
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FIGURE  2-lk.  LOAD  LINE  ANALYSIS 
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MAIN  PERFORMANCE  ANALYSIS  DRIVER  PROGRAM 

COMPUTES  BATTERY  EFFICIENCY 

GENERATES  "INSTANTANEOUS"  l-V  PLOTS 

COMPUTES  ENERGY  STORAGE  GROUP  l-V  CHARACTERISTICS 

COMPUTES  BATTERY  CURRENT  AND  VOLTAGE  CHARACTERISTICS 

COMPUTES  OPERATING  POINT  CURRENT  AND  VOLTAGE 

COMPUTES  POWER  CONDITIONING  AND  DISTRIBUTION  GROUP 
CURRENT  AND  VOLTAGE  CHARACTERISTICS 

COMPUTES  POWER  SOURCE  GROUP  l-V  CHARACTERISTICS 

COMPUTES  SHUNTER  LIMITER  l-V  CHARACTERISTICS 
GENERATES  OUTPUT  PLOTS  AND  TABLES 


FIGURE  2-15.  PERFORMANCE  ANALYSIS  SUBPROGRAM  ELEMENTS 


The  second  form  of  graphical  output  is  shown  in  Figures  2-22  and  2-23. 
This  output  includes  separate  graphical  illustrations  of: 

• Power  Source  Group  Summary  Plots. 

• Battery  Performance  Summary  Plots. 


These  illustrations  show  the  temporal  trends  in  selected  data  taken  from 
the  tabular  presentations  (Figures  2-11  to  2-18)  discussed  above.  In 
particular,  the  information  in  Figure  2-23  shows  that  available  battery 
capacity  (a  critical  performance  parameter)  appears  to  be  holding  steady 
(for  the  design  represented  by  these  curves)  during  the  low  solar  insola- 
tion winter  days. 
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FIGURE  2-16.  UNREGULATED  BUS  SUMMARY  DATA 
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FIGURE  2-18.  ENERGY  STORAGE  UNIT  SUMMARY  DATA 
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FIGURE  2-20.  INSTANTANEOUS  CURRENT-VOLTAGE  PLOTS 
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2-22.  POWER  SOURCE  GROUP  SWIMARY  PLOTS 
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3.  PROGRAM  DATA 


Exact  values  of  all  numerical  program  data  may  be  found  in  Program  Docu- 
mentation (Volume  II,  Users  Manual).  The  basis  for  key  elements  of  this 
data,  however,  will  be  presented  in  this  section  of  the  report. 

Program  Data  may  be  conveniently  placed  in  three  categories  corresponding 

4 

to  the  three  major  categories  of  equipment,  i.e., 

• Power  Source  Group. 

• Energy  Storage  Group. 

• Power  Conditioning  and  Distribution  Group. 

Key  data  for  these  equipment  groups  is  discussed  below. 

3.1  Power  Source  Group 

Information  needed  in  determining  Power  Source  Group  performance  and  char- 
acteristics includes  estimates  of  solar  insolation.  These  estimates  may  be 
obtained  by  the  methods  indicated  in  Appendix  B.  This  appendix  also  con- 
tains the  following  data  needed  in  making  estimates  of  solar  insolation: 

• Time  Zones  in  the  United  States  (Figure  B-2). 

• Time  Zone  Numbers  in  United  States  for  Standard  Time  (Table  B-l). 

• Clearness  numbers  of  Nonindustrial  Atmosphere  in  the  United  States 
(Figure  B-3). 

• Solar  Radiation  Variables  (Table  B-2). 

• Solar  Radiation  Fourier  Coefficients  (Table  B-3). 

• Cloud  Cover  Modifiers  (Table  B-H). 

• Cloud  Cover  Modifier  Equations  (Table  B-5). 

Another  important  parameter  affecting  performance  is  the  solar  cell  spec- 
tral correction  factor.  Data  received  from  the  U.S.  Coast  Guard  (Ref- 
erence 5)  indicate  that  the  short  circuit  current  of  a solar  cell,  above 
the  atmosphere,  is  related  to  solar  cell  short  circuit  current  at  the 
earth's  surface  by  the  following: 


3-1 


50^0-27 


I 

SC 


p 

( 135.6  mw/ cm  I 
I AMO  1 


= (1.11+6)  I 

SC 


2 

100  mw/cm  | 
AMI  i 


(3-1) 


where : 


I 

sc 


1 135*6  mw/cm^  | 
| AMO  i 


Solar  Cell  Short  circuit  current  above  the 
atmosphere  (AMO)  at  a standard  solar  insolation 
of  135.6  milliwatts /cm2 


I 


p 

| 100  mw/cm  I 
sc  I AMI  I 


Solar  Cell  Short  circuit  current  at  the  earth's 
surface  (AMl)  at  a standard  solar  insolation  of 
100  milliwatts/cm2 


Thus,  the  spectral  correction  factor  (SPECOR)  is: 

specok  = H§r)(rlfer)  - ia83 


(3-2) 


Variations  in  the  average  daily  temperature  are  typically  represented  by 
the  yearly  temperature  profile  for  Los  Angeles,  California  (Reference  6), 
shown  in  Figure  3-1. 


Information  on  the  performance  of  terrestrial  solar  arrays  under  test  were 
obtained  from  the  Coast  Guard  (Reference  7).  The  data  obtained  from  these 
tests  were  used  to  estimate  the  equivalent  solar  cell  characteristics  for 
two  types  of  arrays  of  interest  to  the  Coast  Guard.  These  cell  characteris- 
tics are  shown  in  Figures  3-2  and  3-3.  Extrapolation  of  cell  performance 
beyond  measured  data  is  required  because  of  the  DSPA  Computer  Program 
methodology  used  for  solar  array  performance  characteristics. 

A number  of  solid-state  devices  have  been  used  as  solar  array  isolation 
diodes  by  the  Coast  Guard  (Reference  7).  The  characteristics  of  two  of 
these  devices  are  shown  in  Figures  3-^+  and  3-5.  One  of  these  devices  is  a 
rectifier  diode  (Figure  3-M . The  other  device  is  a zener  diode  which 
was  used  by  the  Coast  Guard  as  a rectifier  diode  because  of  its  low  reverse 
leakage  current  (at  less  than  zener  breakdown  voltage  - 39  volts). 
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FIGURE  3-3.  EQUIVALENT  SOLAR  CELL  CHARACTERISTICS 
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FIGURE  3-U.  SOLAR  ARRAY  ISOLATION  DIODE  CHARACTERISTICS 
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FIGURE  3-5.  SOLAR  ARRAY  ISOLATION  DIODE  CHARACTERISTICS 

Shunt  limiters,  included  in  the  DSPA  Computer  program  are: 

• Ordinary  zener  diodes. 

• Temperature  compensated  zener  diodes, 

• Active  shunt  limiter. 

Temperature  coefficients  for  ordinary  zener  diodes  are  shown  in  Figure  3-6 
(Reference  8).  The  information  in  this  figure  shows  that  zener  diodes 
with  a breakdown  voltage  in  the  vicinity  of  5-1  volts  have  a zero  tempera- 
ture coefficient  and  may  also  be  considered  "temperature-compensated." 

The  impedance  of  ordinary  zener  diodes  is  shown  in  Figure  3-7.  The  data 
shown  in  Figures  3-6  and  3-7  are  used  to  estimate  the  current-voltage  char- 
acteristics of  ordinary  zener  diodes. 

Characteristics  of  a typical  series  of  temperature-compensated  zener 
diodes,  as  a function  of  temperature,  are  shown  in  Figure  3-8  (Reference  8). 
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♦25°C 


5040-27 

The  performance  characteristics  of  an  active  shunt  limiter  used  by  the 
Coast  Guard  (Reference  7)  are  given  in  Figure  3-9* 

3.2  Energy  Storage  Group 

Data  have  been  obtained  for  two  storage  cell  types  of  interest  to  the  Coast 
Guard.  The  batteries  of  interest  are: 

• Lead  - Acid  Batteries. 

• Ilickel  - Cadium  Batteries. 

The  effect  of  temperature  and  discharge  rate  on  lead-acid  battery  capacity 
is  shown  in  Figure  3-10  (References  9-15).  As  shown  in  this  figure,  stan- 
dard discharge  conditions  are: 

• Normalized  Discharge  Rate  - 0.05  Ampere s/Anpere-Hours  (20  hour  rate). 

• Temperature  - 20 °C. 
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FIGURE  3-9.  ACTIVE  SHUNT  LIMITER  CHARACTERISTICS 
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FIGURE  3-10.  LEAD  - ACID  BATTERIES  VARIATION  III  DISCHARGE 
CAPACITY  WITH  TEMPERATURE  ADD  DISCHARGE  RATE 


Since  the  normalized  discharge  rate  for  operation  on  a buoy  is  as  low  as 
0.0005  hr  \ the  data  in  Figure  3-10  indicate  that  battery  discharge 
capacities  as  high  as  50  percent  greater  than  the  capacity  under  standard 
discharge  rates  may  be  utilized  for  these  Coast  Guard  applications. 


Current-Voltage  characteristics  of  lead-acid  batteries  at  20°C  are  shown 
in  Figure  3-11.  (References  9,  10,  13,  1^,  15)«  The  voltages  in  Fig- 
ure 3-11  may  be  adjusted  for  other  temperature  by  use  of  the  following 
(Reference  9): 


'V  - V 

C CO 


0.0035  volts/ °F 


NORMALIZED  CHARGE  DISCHARGE  RATE,  houri 
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TEMPERATURE,  20*C 


STORAGE  CELL  VOLTAGE,  volti/CELL 

FIGURE  3-11.  LEAD-ACID  BATTERY,  CURRENT  - VOLTAGE  CHARACTERISTICS 


where 

V = Storage  Cell  Voltage  at  Tc , volts/cell 
c 

V = Storage  Cell  Voltage  at  T^,  volts/cell 
T^  = Storage  Cell  Temperature,  °F 

T = Storage  Cell  Reference  Temperature,  °F 

co 

The  instantaneous  charge  efficiency  of  lead-acid  batteries  is  given  in 
Figure  3-12.  (References  5,  9,  12).  A characteristic  of  these  storage 
cells  is  the  decrease  in  charge  efficiency  with  increasing  charging 
currents . 

The  effect  of  depth-of-discharge  on  lead-acid  battery  cycle  life  is  shown 
in  Figure  3-13.  (Reference  11).  The  data  presently  used  in  the  DSPA 
Computer  Program  is  also  shown.  This  data  is  only  an  engineering  estimate 
and  should  be  updated  as  the  latest  information  from  the  HAD-Crane  tests  on 
lead-acid  batteries  is  received. 
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J-  = NORMALIZED  CHARGE  CURRENT,  -2!!£2£ 

v-  ampere-hour 


FIGURE  3-12.  LEAD  - ACID  BATTERIES,  INSTANTANEOUS  CHARGE  EFFICIEIICY 


FIGURE  3-13-  CHARGE/ DISCHARGE  CYCLE  LIFE  OF  Pb  - ACID  BATTERIES 
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Maximum  allowable  charging  currents  for  lead-acid  Batteries  are  shown  in 
Figure  3-lb  (Reference  15).  The  allowable  current  values  are  used  to 
estimate  the  current  cutoff  on  a battery  charger. 

Electrolyte  characteristics  for  lead-acid  batteries  are  illustrated  in 
Figures  3-15  and  3-16  (Reference  15).  Both  electrolyte  specific  gravity 
and  freezing  point  are  shown.  As  the  data  in  Figure  3-10  shows,  when  the 
electrolyte  temperature  approaches  the  freezing  point,  battery  discharge 
capacity  falls  off  sharply. 

Host  of  the  information  on  nickel-cadmium  (Hi-Cd)  batteries,  supplied  with 
the  DSPA  Computer  Program,  was  based  on  the  results  of  tests  on  Hi— Cd 
storage  cells  used  for  the  Viking  Orbiter  spacecraft  (Reference  3). 

Typical  information  used  in  determining  appropriate  Hi— Cd.  battery  perform- 
ance is  shown  in  Figures  3-17  and  3-l8  (Reference  l6).  Current-voltage 


' 10  100 

IATTERY  OtSCMARGF  OPACITY.  hour*  (At  80* F , 20-hour  RATE) 


- -.  ’V  • " :*/•  • ' LRY  CHARGE  'ilRRKHT,  LEAD  - ACID  BATTERIES 
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CELL  VOLTAGE,  volt. 

FIGURE  3-17.  NICKEL  - CADMIUM  BATTERIES,  CURRENT  - VOLTAGE  CHARACTERISTICS 
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INSTANTANEOUS  EFFICIENCY  OF 
HERMETICALLY  SEALED  NICKEL- 
CADMIUM  CELLS  AS  A FUNCTION 
OF  STATE  OF  CHARGE 

NORMALIZED  CURRENT  AS  A 
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STATE  OF  CHARGE  (Q)  (%  OF  NAMEPLATE  CAPACITY) 

FIGURE  3-18.  NICKEL  - CADMIUM  BATTERIES,  CHARGE  EFFICIENCY 
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characteristics  of  Ni-Cd  storage  cells  are  shown  in  Figure  3-17*  Instan- 
taneous charge  efficiencies  of  these  cells  are  illustrated  in  Figure  3-18. 

The  data  in  Figure  3-18  show  that  the  charge  efficiency  of  Ili— Cd  batteries 
increases  as  charge  rate  increases  (the  opposite  of  lead-acid  batteries). 

Cycle  life  of  Ni-Cd  batteries  for  use  in  design  synthesis  is  given  in  Fig- 
ure 3-19  (Reference  17).  I 

The  characteristics  of  a series  regulator  used  as  a battery  charger  by 
the  Coast  Guard  is  shown  in  Figure  3-20  (Reference  5).  This  regulator  has 
no  current  limit  to  prevent  excessive  charging  currents. 

3.3  Power  Conditioning  and  Distribution  Group 

The  load  requirements  on  the  power  system  are  based  on  various  lamp/ flasher 
combinations.  Two  types  of  lamps  are  under  consideration:  tungsten  fila- 
ment lamps  which  are  currently  in  use  and  xenon  flash  lamps  (over  a cur- 
rent range  which  yields  illumination  levels  equal  to  the  tungsten  lamps). 
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FIGURE  3-19.  CHARGE/DISCHARGE  CYCLE  LIFE  OF  NICKEL  - CADMIUM  BATTERIES 
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AVERAGE  NAD-CRANE  DATA 
1.5  hr  ORBIT  BASED  ON  TEST 
TO  FAILURE  INFORMATION 

SAFT  8 hr  ORBIT  20"C 

HUGHES  TEST  ON  GE  15  ohCELLS 
8 hr  ORBIT.  NO  FAILURES  AFTER 
451  CYCLES.  TESTING  STARTED 
1976 

NAD-CRANE  TEST  1.5  hr  ORBIT 
QUL  1 ON  9.8  ML  1964 


1.5  hr  CYCLE 
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REGULATOR  INPUT  VOLTAGE,  vdc 

FIGURE  3-20.  BATTERY  CHARGER  (REGULATOR)  CHARACTERISTICS 

Flasher  types  are  those  presently  under  development  or  operational  with  the 
Coast  Guard  (Reference  18),  i.e.: 

• Type  I - Simple  Flashing. 

• Type  II  - Group  Flashing. 

• Type  III  - Complex  Flashing. 

• Type  IV  - Fixed. 

These  flashers  have  the  characteristics  indicated  below: 


'Type  I - Simple  flashing  consisting  of  a flash  of  F seconds, 
with  an  eclipse  of  seconds,  the  total  period  being  + F^ 


alternating 
seconds . 


The  most  common  Type  I characteristics  are  specified  in  Table  3-1. 


- Group  flashing  consisting  of  a group  of  II  flashes  all  of  the  same 


duration  F , with  the  inter-flash  eclipse  of  duration  E^,  the  intergroup 
eclipse  of  duration  E , and  the  overall  period  of  II F ^ + (II  - l)  E^  + E^. 
Common  Type  II  characteristics  are  in  Table  3-2. 
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TABLE  3-1. -TYPE  I FLASHER  CHARACTERISTICS 


F1 

E1 

Characteristic  Hrkg. 

Duty  Cycle 

0.3  seconds 

0.7  seconds 

QK  FL 

(0.5) 

0.300 

0.5 

2.0 

FL  2.5 

(0.5) 

0.200 

0.4 

3.6 

FL  4 

(0.4) 

0.100 

1.0 

5.0 

FL  6 

(1.0) 

0.167 

3.0 

3.0 

E HIT  6 

(3.0) 

0.500 

3.0 

1.0 

0CC  4 

(3.0) 

0.750 

TABLE  3-2. -TYPE  II  FLASHER  CHARACTERISTICS 


II 

F1 

Ei 

E2 

Characteristic 

Marking 

6 

0.3  seconds 

0.7  seconds 

4.7  seconds 

I QK  FL  (6X0.3) 

0.180 

2 

0.U  seconds 

0.6  seconds 

3.6  seconds 

GP  FL  5 (2X0.4) 

0.160 

TYPE  III  - Complex  flashing  consisting  of  any  other  combination  of  short 
and  long  flashes  ( F^ , F^,  etc.)  and  short  and  long  eclipses  (E  , E0,  etc.). 
The  only  common  characteristic  of  this  type  is  the  Morse  Code  "A,"  indi- 
cated in  Table  3-3. 


TABLE  3-3. -TYPE  III  FLASHER  CHARACTERISTICS 
(MORSE  CODE  "A") 


F1 

E1 

F2 

E2 

Characteristic 

Marking 

— 

2.0  seconds 

5.0  seconds 

MO (A)(0.4,2,0) 

0.300 

L8 


r : i 


Type  IV  - A flasher  with  a fixed  characteristic  will  provide  all  the 
operational  features  of  the  other  types  but  will  allow  the  lamp  to  burn 
steadily  with  no  time-coding  imposed  on  the  voltage.  The  characteristic 
marking  is  "FIX"  and  the  "duty  cycle"  marking  is  "1.000." 

The  data  shown  above  have  been  incorporated  in  the  DSPA  Computer  Program. 
The  program  user  may  then  simply  select  a desired  flasher  by  using  the 
code  number  indicated  in  Table  3-M 

Nominal  currents  (for  Tungsten  Lamps)  are  in  the  range  0.25  amperes  - 
3.05  amperes  (Reference  18).  The  nominal  current  must  be  adjusted  for 
cold  filament  surge  at  lamp  turn-on.  Values  of  the  cold  filament  surge 
coefficient  are  shown  in  Figure  3-21  (Reference  19). 

Voltage  regulation  characteristics  of  the  flasher  regulator  are  shown  in 
Figure  3-22  (Reference  l8).  This  regulator  is  only  under  load  when  the 
lamp  is  turned  on. 


TABLE  3-M -FLASHER  CODE  NUMBERS 


IFTYPE  (CODE  NUMBER) 

CHARACTERISTIC  MARKING 

1 

QKFL 

(0.3) 

2 

FL  2.5 

(0.5) 

3 

FL  L 

(o.M 

1* 

FL  6 

(1.0) 

5 

EIIIT  6 

(3.0) 

6 

OCC  U 

(3.0) 

7 

IQKFL 

(6  x 0.3) 

8 

GPFLS 

(2  x 0.U) 

9 

MO  ( A ) 

(0.U.2.0) 

10 

FIX 
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Characteristics  of  the  "housekeeping"  load  regulator  are  shown  in  Fig- 
ure 3-23  (Reference  5)-  These  "housekeeping"  loads  are  always  effective 
and  represent  a permanent  load  on  the  power  system. 


FIGURE  3-23.  HOUSEKEEPING  LOADS  - AIDS  TO  NAVIGATION 
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I1+.  PROGRAM  VERIFICATION 

The  final  phase  in  the  development  of  the  DSPA  Computer  Program  was 
program  verification.  For  this  phase,  the  Coast  Guard  supplied  JPL  with 
performance  data  on  five  actual  solar  array/battery  power  system  designs. 

| 

The  purpose  of  program  verification  was  to  provide  a check  on  the  valid- 
ity of  the  performance  analysis  segment  of  the  program  by  comparison  of 
the  predicted  operational  characteristics  of  a known  power  system  with  the 
actual  operational  characteristics  of  that  same  power  system.  The  power 
systems  against  which  the  comparisons  were  made  were  selected  from  those 
presently  undergoing  instrumented  tests  at  the  Coast  Guard  R&D  Center. 
Recorded  values  of  all  critical  parameters,  from  these  tests,  were  made 
available  to  JPL  by  the  U.S.  Coast  Guard  (References  5»  12). 

Major  characteristics  of  the  five  power  systems  used  in  this  comparison 
are  shown  in  Table  1+-1;  while,  the  general  arrangement  of  the  power 
system  equipment  is  shown  in  Figure  4-1.  The  information  in  Table  4-1 
indicates  that,  to  all  intents  and  purposes.  Systems  No.  2 6 and  29  are 
identical . 

Another  key  performance  parameter  measured  by  the  Coast  Guard  was  solar 
array  current  degradation  — shown  in  Figures  4-2  to  4-6.  Relatively 
erratic  performance  was  demonstrated  by  the  solar  arrays  on  Systems  2 and 
9.  The  System  47  solar  array  also  experienced  erratic  performance  due  to 
corrosion  effects  early  in  the  test.  In  fact,  the  only  two  arrays 
declared  successful  by  the  Coast  Guard  for  purposes  of  verification  were 
those  used  for  Systems  26  and  29. 

The  major  parameter  used  to  determine  successful  prediction  of  power 
system  operation  was  battery  state-of-charge . The  reason  for  using 
this  parameter  was  that  if  battery  charge  drops  to  zero  at  any  time 
during  a mission  the  power  system  is  considered  to  have  failed. 
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TABLE  lt-l.-NAVAID  POWER  SYSTEM  CHARACTERISTICS 


System 

Number 

Solar  Array 

Array 

Isolation 

Diode 

Shunt 

Limiter 

Battery 

Lamp 

(amp) 

O 

Heliotek 

HGL-V-3785-03 

1cm  x 2cm  P on  N cells* 

1N131*1A 

Heliotek 

Globe; 

12  volts 

20  amp-hours/ 
battery 

0.55 

1*0  cells-in-series 

15  cells-in-parallel 

3 batteries- 
in-parallel 

9 

Heliotek 

1N13U1A 

Heliotek 

Gates ; 

2.1  volts/cell 

0.55 

HGL-V-3785-03* 

6 cells-in- 
series 

Same  as  system  No.  2 

5 amp-hours/ 
battery 

12  batteries- 
in-parallel 

2 6 

Heliotek 

HW-3785-Ol* 

1cm  x 2cm  P on  N cells* 

33  cells-in-series 

15  cells-in-parallel 

1N13U1A 

WISCO; 

12  volts 
(2  batteries- 
in-series ) 

100  amp-hours/ 
battery 

0.77 

29 

Heliotek 

HW-3785-02* 

Same  as  system  No.  26 

1N131*1A 

Same  as  system 
No.  26 

0.77 

'See  Figure  3-2. 
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TABLE  1+-1.-NAVAID  POWER  SYSTEM  CHARACTERISTICS  (contd) 


System 

Number 


olar  Array 


Centralab  No.  17 

2cm  x 2cm  P on  N 
cells** 

1+8  cells-in-series 
6 cells-in-parallel 


Array 

Isolation  Shunt 
Diode  Limiter 


1N2992 


Battery 

Lamp 

(amp) 

Globe; 

12  volts 

20  amp-hours/ 
battery 

2 batteries- 
in-parallel 

0.55 

*See  Figure  3-3 
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FIGURE  1+-1.  POWER  SYSTEMS  TESTED  BY  U.S.  COAST  GUARD 
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FIGURE  h-2.  SOLAR  ARRAY  CURRENT  DEGRADATION  SYSTEM  NO.  2 
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FIGURE  b-b.  SOLAR  ARRAY  CURRENT  DEGRADATION  SYSTEM  NO.  2 6 


FIGURE  l*-5.  SOLAR  ARRAY  CURRENT  DEGRADATION  SYSTEM  NO.  29 
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FIGURE  U-6.  SOLAR  ARRAY  CURRENT  DEGRADATION  SYSTEM  HO.  U7 


The  results  of  the  verification  tests  of  the  DSPA  Computer  Program  are 
shown  in  Figures  U-7  to  U-ll.  Each  of  these  figures  shows  two  values  for 
the  battery  storage  capacity.  The  standard  capacity  (CBSTD)  was  taken 
from  Table  U-l  (except  in  the  case  of  System  9 where,  for  convenience, 
the  battery  capacities  were  doubled,  and  the  batteries-in-parallel  were 
halved).  These  values  represent  battery  capacity  at  the  standard  dis- 
charge rate  of  0.05  ampere/ampere-hour  (see  Figure  3-10).  However,  as 
the  battery  discharge  rate  decreases,  the  storage  capacity  of  the  battery 
increases.  Hence,  the  "actual"  battery  capacities  (CB)  are  also  shown 
in  Figures  U-7  to  b-11.  These  values  of  CB  represent  the  battery 
capacities  used  as  input  to  the  DSPA  computer  program  which  enabled  the 
program  to  predict  the  battery  state-of-charge  (as  a function  of  time) 
shown  in  Figures  U-7  to  U-ll. 

The  predicted  battery  states-of-charge  for  Systems  2 and’  9 were  still 
a fair  approximation  of  the  measured  state-of-charge  in  spite  of  the 
erratic  performance  of  tne  solar  arrays  for  these  systems. 
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FIGURE  4-9.  COMPUTER  PROGRAM  VERIFICATION  TEST  SYSTEM  NO.  26 


COMPUTER  PROGRAM  VERIFICATION  TEST  SYSTEM  NO.  29 
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FIGURE  4-10 . 


STATE-OF-CHARGE 
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FIGURE  1+-11.  COMPUTER  PROGRAM  VERIFICATION  TEST  SYSTEM  NO.  1*7 

Predicted  battery  states-of-charge  for  Systems  2 6 and  29  were  very  close 
to  the  measured  states-of-charge.  These  predictions  made  with  success- 
ful power  systems  are  a true  measure  of  the  capability  of  the  DSPA  com- 
puter program  to  simulate  power  system  operation  over  a one  year  period. 

The  final  verification  test  performed  with  System  1+7  data  only  shows  a 
general  trend  agreement  between  measured  and  predicted  battery  state-of- 
charge.  This  is  due  to  the  initial  corrosion  failure  of  the  solar  array 
with  the  attendant  erratic  operation.  This  comparison  is  not  consid- 
ered a true  test  of  the  DSPA  Computer  Program  capabilities. 
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APPENDIX  A 

BALANCED  POWER  SYSTEMS 


LIST  OF  FIGURES 
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A schematic  of  the  Navigation  Aid  Power  system  is  shown  in  Figure  A-l. 

There  are  three  operational  modes  for  this  power  system,  i.e.: 

• Power  Source  Group  Only. 

• Battery  Only. 

• Share-Mode . 

When  there  is  no  power  available  from  the  Power  Source  Group  (PSG),  then 
the  loads  are  suj plied  only  by  the  battery.  Hence,  the  energy  balance  is: 

^PPCD^J  ^At^J  = nD  ^PBD^J  lA"l) 

where : 

(At)T  = Incremental  load  periods  during  Power  Source  Group  shutdown, 
hours 

(P  ) = Power  Conditioning  and  Distribution  Group  (PCDG)  loads  during 

u ^ PSG  shutdown,  watts 

(P  )T  = Battery  Discharge  Power  during  PSG  shutdown,  watts 
BD  J 

nD  = Efficiency  of  battery  discharge  device 
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During  Share-Mode  operation  when  the  loads  are  supplied  by  both  the  Power 
Source  Group  and  the  battery,  the  energy  balance  is: 

^PPCD^K  = PPSG  + nD  ^PBdV  (A-2; 

where : 


(At)^  = Incremental  load  periods  when  PSG  is  charging  the  battery  and 
powering  the  loads,  hours 

n„TI  = Battery  charger  efficiency 

On 

(PDm,)T  = Battery  Charge  Power  during  charging  periods,  watts 

13  On  L 

(P  ) = PCDG  loads  during  battery  charge  periods,  watts 

r OJJ 

During  the  PSG  shutdown  periods: 


5040-27 


During  PSG  and  Battery  Share-Mode  periods : 


^BD^fC  “ ^PBD^K  ^At^K 


(A-5 ) 


where:  (E^),.  = Battery  energy  discharged  during  Share-Mode  periods, 

dU  K-  . , , 

watt-hours 


Combining  Eqs.  (A-i)  and  ( A- 4 ) : 


^BD^J  " 


(Wj  (At). 


(A -6) 


Combining  Eqs.  (A-2)  and  (A-5): 


^BD^K  = 


^PPCdV  (AtV  PPSG  tAt‘>K 


( A— 7 ) 


The  total  battery  discharge  energy  is: 


r'BDT  ^EBD^J  + (EBD^K 


( A— 8 ) 


where:  EDnm  = Total  battery  energy  discharged  during  "mission,"  watt-hours 

DJJl 


Combining  Eqs.  (A-6),  ( A— 7 ) and  ( A— 8 ) : 


'' ' PCp)  J ^At^J  ^ PPCD^K  (At^K  PPSG  ^AtV 
EbDT  nD  "o'  "d 


During  period  of  battery  charge: 


(A-9) 


EBCHT  = ^PBCH^L  (At)L 


where:  = Total  battery  charge  energy  during  "mission,"  watt-hours 

BCni 


(a-io: 
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Combining  Eqs.  (A-3)  and  (A- 10): 


SBCHT  ^nCH^  ^PPSG^  ^At^L  “ ^nCH^  ^PPCD^L  ^At^L 


(A-ll) 


A balanced  power  system  is  defined  as  one  wherein  the  PSG  supplies  just 
enough  energy  to  the  Energy  Storage  Group  (ESG)  to  replace  that  energy  dis- 
charged by  the  ESG.  Thus: 


EBCHT  ~ VBE  EBDT 


(A-12) 


where:  = Battery  energy  charge/discharge  ratio  (inverse  of  battery  energy 

charge  efficiency) 

Combining  Eqs.  (A-9),  (A-ll)  and  (A-12)  and  rearranging,  the  average  PSG 
power  level  requirement  is  then: 

_ fvBE5  f(PPCD)J  + (PPCD^K  + fnDnCH5  EPPcdL  (At  1 

PPSC <VCH>  (4tlL  * <V>  ' 3) 


Defining: 


^PCD^J  ^PPCD^J  ^At^J 


^epcdV  " ^ppcdV  ^AtV 


^EPCD^L  " ^ PPCD^K  ^At^L 


Load  energy  ( from  PCD  group ) 
during  each  of  the  applicable 
power  system  operational 
modes,  watt-hours 


(A-lU) 


(A-15) 


(A- 16) 


Combining  Eqs.  (A-13),  (A-l^),  (A-lp)  and  ( A— l6 ) : 


^VBE^  jEPCP^J  * ^EPCdV.  * ^nDnCH^  ^PCD^L 
^nDnCH^  ^At^L  + ' VBE ^ ^AtV 


(A-17) 
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The  total  energy  required  of  the  PSG  is  thus: 


EPSG  ^PPSG^  + ^At^L 


(A-18) 


where:  EpgG  = Power  Source  Group  Total  energy  requirement  for  the  "mission,' 


watt-hours 


A-7/A-8 
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"Procedure  for  Determining  Heating  and  Cooling  Loads  for  Computerized  Energy 
Calculations"  - American  Society  of  Heating,  Refrigerating  and  Air  Condition- 
ing Engineers;  1971. 


B-l/B-2 
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SUIJ 


An  Algorithm  to  Find  Solar  Position,  and  Intensity  of  Solar 
Radiation  Incident  on  the  Outer  Surfaces  of  Buildings 


FIGURE  B-l.  DEFINITIONS  OF  ANGLES 

INPUT: 


L 


l 


TZN 


d 

t 


g 

CCM 


Latitude,  degrees, 


Longitude,  degrees, 


I- +North*| 
|_»Soutb  J 

[+V.’est*J 
-East J 


Tine  2one  number  (hours  behind  Greenwich  mean  time)  (see 
Figure  B-2  and  Table  B-l) 

Date,  days  (from  start  of  year),  (l  - 3 66) 

Time,  hours  (after  midnight),  (0  - 2M 

Ground  reflectivity 


Cloud  cover  modifier  (Output  of  CCM) 


B-3 


— 


Mi  LUMBERS  III  ThL  U.C.  F R 


rface  azimuth  angle,  degrees  (from  South) 


Depending  upon  type  of  industry  in  the  locality,  use  Clearness  Number  between 
O.T  and  0.9.  Otherwise,  use  the  Clearness  Numbers  of  nonindustrial  atmosphere 
which  are  given  in  Figure  B-3. 


TIME  ZONE 

TT 

1 Cut 

Atlantic 

J* 

F astern 

5 

'entral 

6 

tl  >untain 

T 

! acific 

8 

S SUMME  RlMjJASn 
A *iNTf  N IN  l)  J » MA| 


F JURE.  -3.  'LEARNESS  NUMBER;'  OF  II 
A1:  M .TH!  HE  IN  UNITED  STATE:' 


mri  • and  sunset  time,  I:  *urs  (after  midnight) 


ines  of  direct  f.  lar  bear 


ar  radiation 


uae 


lar  azimuth,  degrees 


diffuse  sky  radiation  on  horizontal  surface) 


B :'ky  brightness 
Btu  per  (hr)(s 


BG  Ground  brightness  (=  diffuse  ground  reflected  radiation 
Btu  per  (hr)(sq  ft) 


IDN  Intensity  of  direct  normal  solar  radiation.  But  per  (hr) 
(sq  ft) 
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I Intensity  of  total  solar  radiation  incident  on  surface  (wall), 
Btu  per  (hr) (sq  f t ) 


d,  sky- 


intensity  of  sky  diffuse  radiation  incident  on  surface  (wall), 
Btu  per  (hr)(sq  ft) 


d, ground 


Intensity  of  ground  diffuse  radiation  incident  on  surface 
(wall),  Btu  per  (hr)(sq  ft) 


Table  B-2  lists,  as  function  of  date  and  for  the  northern  hemisphere  only, 
five  variables  related  to  solar  radiation.  These  variables  are  declination 
angle  6;  the  equation  of  time,  ET;  the  apparent  solar  constnat,  A;  the  atmos- 
pheric coefficient,  B;  and  sky  diffuse  factor,  C. 


Table  B-2  could  be  stored  in  the  computer  memory,  but  this  would  necessi- 
tate an  interpolation  procedure.  In  order  to  avoid  such  a problem  and  to  save 
computer  core  space,  6,  ET,  A,  B and  C are  expressed  in  Fourier  Series  form, 
and  the  values  are  calculated  as  a function  of  the  day  of  the  year,  d,  from 
the  following  truncated  Fourier  series. 

* w * d)  + A,  *Cos(3*w*d) 

* w * d)  + B.  * Sin( 3 * u * d) 


) 


ET;  -0  ■ Ax  * Cos  (cj  * d)  + A2  * Cos  (2 
A 

Bj  + E1  * sin(w  * d)  + B2  * Sin(2 


where  w = 2 * tt/366. 


TABLE  B-2. -VALUES  OF  6,  ET,  A,  B AND  C 


Date 

6 

Degrees 

ET 

Hours 

A 

Btu  per 
(hr ) ( sq  ft) 

B 

Air  Mass 

C 

Jan.  21 

-20.0 

-0.190 

390 

0.ll*2 

0.058 

Feb.  21 

-10.8 

-0.230 

385 

0.  lUl* 

0.060 

Mar.  21 

0.0 

-0.123 

376 

0.156 

0.071 

Apr.  21 

11.6 

0.020 

360 

0.180 

0.097 

May  21 

20.0 

0.060 

350 

0.196 

0.121 

June  21 

23.1*5 

-0.025 

31*5 

0.205 

0.13k 

July  21 

20.6 

-0.103 

3l*l* 

0.207 

0.136 

Aug.  21 

12.3 

-0.051 

351 

0.201 

0.122 

Sept.  21 

0.0 

0.113 

365 

0.177 

0.092 

Oct.  21 

-10.5 

0.255 

378 

0.160 

0.073 

Nov.  21 

-19.8 

0.235 

387 

0. 11+9 

0.063 

Dec.  21 

-23.1*5 

0.033 

391 

0.11*2 

0.057 

B-6 
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The  proper  Fourier  coefficients  are  given  in  Table  B-3 


TABLE  B-3. -FOURIER  COEFFICIENTS 


0.302 

0.0 

368.UL 

0.1717 

0.0905 


-0.229 
-0.05 
-l.ll* 
0.0032 
O.OUlO  I 0.0073 


0.0015 


CALCULATION  SEQUENCE: 


1.  Determine  or  calculate  6,  ET,  A,  B,  C and  CCM. 

2.  Calculate  sunrise  and  sunset  time,  h',  radians. 

h'  = Cos  1(-Tan(L)  * Tan(6)) 

Y = h'  * (12/tt) 

Sunrise  time,  SRT,  hr 

SRT  = 12  - Y - ET  - TZN  + SI/ 15 
Sunset  time,  SST,  hr 

SST  = 2l+  - SRT 

3.  Calculate  hour  angle,  h,  degrees. 

h = 15  * (t  - 12  + TZN  + ET)  - £ 


B1 

B2 

B3 

3.851 

0.002 

-0.055 

-0.122 

-0.156 

-0.005 

0.58 

-0.18 

0.28 

-0.001*3 

0.0 

-0.0008 

-0.0031+ 

o.oooi* 

-0.0006 

B-7 
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4.  Check  if  |h | > |h* | . 


If  yes,  make  1=0  and  skip  following  steps. 

If  no,  go  to  next  step. 

5.  Calculate  direction  cosines  of  direct  solar  beams. 

Cos(Z)  = Sin(L)  * Sin(6)  + Cos(L)  * Cos(6)  * Cos(h) 

Cos(W)  = Cos(6)  * Sin(h) 

Cos ( S ) = (1  - (Cos(Z))**2  - (Cos(W) )*#2)**0.5 

If  Cos(h)  > (Tan(6 )/Tan(L) ) , Cos(S)  is  positive,  otherwise  it  is 
negative. 

6.  Calculate  solar  altitude  and  solar  azimuth. 

SALT  = Sin_1(Cos(Z) ) 

If  Cos(S)  > 0,  SAZM  = Sin” X( Cos (W)/Cos (SALT) ) . 

If  Cos (S ) < 0,  SAZM  = 180  - Sin-1 ( Cos (W)/ Cos ( SALT) ) . 

7.  Calculate  intensity  of  direct  normal  solar  radiation. 

IDN  = A * CN  * CCM  * Exp(-B/ Cos (Z ) ) 

8.  Calculate  sky  brightness. 

BS  = C * IDN/(CN**2) 

9.  Calculate  ground  brightness. 

BG  = p * (BS  + IND  * Cos(Z)) 
g 


B-8 
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10.  Calculate  direction  cosines  of  normal  to  surface  (wall). 

a = Cos(WT) 

0 = Sin(WA)  * Sin(WT) 
y = Cos(WA)  * Sin(WT) 

11.  Calculate  cosine  of  angle  of  incidence. 

Cos(n)  = a * Cos(Z)  + 0 * Cos(W)  + y * Cos(S) 

12.  Calculate  intensity  of  direct  solar  radiation  incident  on  surface  (wall). 

ID  = IDN  * Cos(n)  If  Cos(n)  > 0,  otherwise  ID  = 0 

13.  Calculate  intensity  of  sky  diffuse  radiation  incident  on  surface  (wall). 

For  Horizontal  Surface: 


Id,sky 


= C * IDN 


For  Vertical  Surface: 

If  Cos(n)  > -0.2 

Y = 0.55  + 0.1*37  * Cos(n)  + 0.313  * (Cos(n)**2) 

Otherwise  Y = 0.1*5 

I.  , = IDN  * (C  * Y + (p  * (C  + Cos ( Z ) ) )/2) 

d,sky  g 

ll*.  Calculate  intensity  of  ground  diffuse  radiation  incident  on  surface  (wall). 


I 


d, ground 


= BG  * ( (l-a)/2) 


B-9 
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15.  Calculate  intensity  of  total  solar  radiation  incident  on  surface  (wall). 

I = (ID  + I,  . + I,  J 

d,sky  d, ground 


CCM* 


Cloud  Cover  Modifier 


INPUT: 


SALT  Solar  altitude  angle,  degrees  (Output  of  SUN) 

CT  Cloud  type 


( 0 cirrus  or  cirrostratus** 
CT  = | 

( 1 stratus 


TC  Total  cloud  amount  at  time  t 


OUTPUT : 


CCM  Cloud  cover  modifier 

Table  B-4^"  lists  the  value  of  cloud  cover  modifier,  CCM,  for  each  combina- 
tion of  cloud  type,  total  cloud  amount,  and  solar  altitude  angle.  Each  column 
of  Table  B-4,  corresponding  to  a given  cloud  type  and  a solar  altitude,  has 
been  expressed  as  a cubic  relationship  between  total  cloud  amount  and  CCM.  The 
resulting  six  equations  are  given  in  Table  B-5. 


*This  subroutine  is  a stop-gap  until  better  information  is  available  for  the 
effect  of  cloud  cover  on  diffuse  sky  radiation. 

**For  any  other  type  cloud,  use  the  mean  of  CCM  for  CT  = 0 and  1. 

TDerived  from  Boeing  Company  Report,  Summary  of  Solar  Radiation  Observation 
D2-90577-l/Dec.  1964. 


L 
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TABLE  B-4. -CLOUD  COVER  MODIFIER 


Cloud  Type 

Cirrus 

Total 

Stratus 

or  Cirrostratus 

Other 

Cloud 

Amount 

SALT  £ 45° 

SALT  > 45° 

SALT  £ 45° 

SALT  >45° 

SALT  £ 45° 

SALT  > 45° 

1 

o.6o 

0.88 

0.84 

1.0 

0.72 

0.94 

2 

0. 60 

0.88 

0.83 

1.0 

0.715 

0.94 

3 

0.58 

0.88 

0.83 

1.0 

0.705 

0.94 

4 

0.58 

0.87 

0.82 

1.0 

0.70 

0.935 

5 

0.57 

0.85 

0. 80 

0.99 

0.685 

0.92 

6 

0.53 

0.83 

0.77 

0.98 

0.65 

0.905 

7 

0.49 

0.79 

0.74 

0.95 

0.615 

0.87 

8 

0.43 

0.73 

0.67 

0.90 

0.55 

0.815 

9 

0.35 

0.6l 

0.60 

0.84 

0.475 

0.725 

10 

0.27 

0.46 

0.49 

1 

0.74 

0.38 

0.42 

TABLE  B- 5. -CLOUD  COVER  MODIFIER  EQUATIONS 


CT  = 

0,  SALT  £ 45° 

CCM  = 0.598  + 0.00026 

* TC 

+ 0.0021  * (TC#»2)  - 

0.00035  * 

(TC**3) 

CT  = 

0,  SALT  > 45° 

CCM  = 0.908  - 0.03214 

* TC 

+ 0.0102  * (TC*#2)  - 

0.00114  * 

(TC**3) 

CT  = 

1,  SALT  £ 45° 

CCM  = 0.849  - 0.01277 

* TC 

+ 0.0036  * (TC«#2)  - 

0.00059  * 

(TC**3) 

CT  = 

1,  SALT  > 45° 

CCM  = 1.010  - 0.01394 

* TC 

+ 0.00553  * (TC»*2) 

- 0.00068  * 

(TC**3) 

CT  # 

0 or  1,  SALT  £ 45° 

CCM  = 0.724  - 0.00652 

* TC 

+ 0.00191  * (TC#*2) 

- 0.00047  * 

(TC**3) 
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ELECTRICAL  CHARACTERISTICS 

LIST  OF  FIGURES 

Figure  Title 

C— 1 Power  Conditioning  and  Distribution  Group  

C-2  Typical  Flasher  Characteristics  

C-3  Lamp-Flasher  Regulator  Characteristics  

C-l  Saturated  Region  Operation  

C-5  Saturated  Pass  Transistor  Characteristics  

C-6  Active  Region  Operation  

C-7  Active  Region  Regulator  Characteristics  

C-8  Lamp-Flasher  Regulator  Input  Characteristics  . . . 

C-9  Housekeeping  Regulator  Input  Characteristics  . . . 

C-10  Lamp-Off  Condition  

C-ll  Lamp-On  Condition  

C-12  Lamp-Flashing  Condition  

C-13  Power  Conditioning  and  Distribution  Group  Input 

Characteristics  


Page 

C-3 

C-3 

C-6 

C-7 

C-7 

C-8 

C-9 

C-10 

C-ll 

C-12 

C-13 

C-13 

C-li 


C-l/C-2 


The  Power  Conditioning  and  Distribution  Group  (PCDG)  consists  of  a Lamp 
Flasher  Regulator  in  parallel  with  a Housekeeping  Regulator  as  well  as  ancillary 
equipment.  A schematic  of  the  PCDG  is  shown  in  Figure  C-l  below: 


FIGURE  C-l.  POWER  CONDITIONING  AND  DISTRIBUTION  GROUP 


Typical  characteristics  of  the  flasher  which  controls  the  lamp  duty  cycle 
are  shown  below  in  Figure  C-2. 
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Based  on  the  characteristics  shown  in  Figure  C-2,  the  lamp  duty  cycle  is: 


nlon 

nlon 

I tL0N(i) 

S tL0II^l) 

i=l 

i=l 

nlon 

N 

L0FF 

tl 

(C-l) 


X tL0t/1^  + ^ tL0FF^J) 


where : 


D = Lamp  duty  cycle 
Li 

t^Qij(i)  = Duration  of  ith  lamp  flash  during  repetitive  period,  , seconds 

t (J)  = Duration  of  Jth  lamp  cutoff  during  repetitive  period,  xT , 

LU  r r • u 

seconds 

Tt  = Repetitive  period  of  flasher,  seconds 
N^Ou  = Total  number  of  lamp  flashes  during 
^LOFF  = nujnber  of  lamp  cutoffs  during 


Correcting  for  cold  filament  surge  at  turn-on: 

|CLg  > 1.0  when:  (0.0  < DL  < 1.0)J 
Dt  - o Tt_  < / 


L LS  LR 


|cLS  = 1.0  when:  DL  = 1.0 


where : 


IT  = Actual  lamp  current,  amperes 
Li 

CTC,  = Filament  surge  correction  factor 
Lb 

IT_  = Rated  lamp  current,  amperes 
Ln 


(C-2) 


Finally,  the  average  lamp  current  is: 


rL  = XL  DL 

where:  I = Average  lamp  current,  amperes 

L 


(C-3) 


The  resistance  rating  of  a lamp  is  given  by: 

R - 

LR  IT„ 


(C-L) 
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where: 

V = Voltage  rating,  VDC 
Ln 

ITD  = Current  rating,  amperes 
LK 

RTT  = Resistance  rating,  ohms 
Ln 

Using  the  actual  lamp  current  (including  turn-on  surge),  the  actual  lamp 
resistance  is: 


rl  = T 


(C-5) 


where : 


I = Actual  lamp  current,  amperes 
L 

R = Actual  lamp  resistance,  ohms 
L 


The  effective  lamp  resistance  (due  to  an  on-off  duty  cycle)  is: 


i 


where : 

I = Average  lamp  current , amperes 
L 

R = Effective  lamp  resistance,  ohms 
L 

The  lamp- flasher  regulator  is  a series  regulator  and  hence  has  the  follow- 
ing characteristics. 


JRI  = IL 


(C-7) 


where:  IDT  = Actual  regulator  input  current,  amperes 

n I 


IRI  = IL 


( C-8) 


where:  IDT  = Average  regulator  input  current,  amperes 

nl 


While  the  voltages  are: 


VL  • F fVRI-  TR> 


I 

L 


( C-9 ) 
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where: 


V = Regulator  input  voltage,  VDC 
ni 

V = Lamp  voltage,  VDC 

Li 

T = Regulator  temperature,  °F 


The  regulator  characteristics  of  the  lamp-flasher  are  shown  in  Figure  C-3. 


ACTIVE  REGION 


SATURATED  REGION 


REGULATOR 
TEMPERATURE  ( 


VRI,  REGULATOR  INPUT  VOLTAGE,  vdr. 


FIGURE  C-3.  LAMP-FLASHER  REGULATOR  CHARACTERISTICS 

When  the  lamp-flasher  regulator  is  operating  in  the  saturated  region 

< VRjgA),  the  regulator  and  lamp  combination  may  be  represented  by 
the  schematic  shown  in  Figure  C-L. 


The  characteristics  of  a saturated  pass  transistor  are  shown  in  Figure  C-5. 


The  voltage  drop  across  the  regulator  is  given  by: 


AVRS  = VRI0  + TL  ZRS 


(c-io: 


c-6 


FIGURE  C-U.  SATURATED  REGION  OPERATION 


AVrs,  VOLTAGE  DROP  ACROSS  "SATURATED"  PASS  TRANSISTOR,  vdc 


FIGURE  C-5.  SATURATED  PASS  TRANSISTOR  CHARACTERISTICS 
But,  by  definition: 

"ns  * vri  ' vl  (c-u) 

Combining  Eqs.  (C-10)  and  (C-ll),  the  regulator  output  voltage  (i.e.,  the 
lamp  voltage)  is  given  by: 


VL  = VRI  - VRIO 


IL  ZRS 


(C-12) 


C-7 
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The  voltage  drop  across  the  lamp  is: 


VL  = JL  RL 


Combining  Eqs.  (C-12)  and  (C-13),  the  lamp  current  is: 


V - V 
RI  RIO 


RL  + ZRS 


(C- 


Substituting  Eq..  (C-7)  into  (C-lU)  yields,  for  the  actual  regulator 
input  current: 


RI 


- /VRI  ~ VRIo\ 

= \ rl+  zrs  ) ; 


(V  < V < V ) 
v RIO  RI  - RISA; 


(C- 


and  by  using  the  same  methods,  the  effective  regulator  input  current  is: 


RI 


/V  ; (V  < 

l »L  + ZES  ) M° 


V < V ) 
RI  - RISA; 


(C- 


where:  I = Effective  regulator  input  current,  amperes 

HI 


When  the  lamp-flasher  regulator  is  operating  in  the  active  region 
(VRI  > VRISA),  the  regulator  and  lamp  combination  may  be  represented  by  the 
schematic  shown  in  Figure  C-6.  ; 
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Voltage  regulator  output  characteristics  in  the  "pass  transistor  active" 
region  are  shown  in  Figure  C-J. 

In  the  active  region,  it  may  be  shown  that  the  regulator  characteristics 
are  given  by: 


LB 


tl  zra 


(C-17) 


where:  = Regulator  output  voltage  at  zero  current,  VDC 


Combining  Eqs.  ( C— 13 ) and  (C-17),  the  lamp  current  is: 


LB 


L RL  + ZRA 


But : 


VLB  = VLB  {VRI}  i-e':  VLB  is  a function  of  Vri 


(C-l8) 


(C-19) 


CURRENT, 

cm  pores 


VALUES  OF  V, 

LAMP  VOLTAGE  (REGULATOR  OUTPUT  VOLTAGE) 


LB 

voc 
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Hence,  combining  Eqs.  (C-7),  (C-l8)  and  (C-19): 


RI 


/vyv 

\ BL  + ZRA  / ’ 


(V  > V ) 
' RI  RISA' 


( C-20 ) 


and  by  using  the  same  methods,  the  effective  currents  are: 


(C-21) 


Finally,  it  should  be  noted  when 


JL  = TRI  = °*°  ) 

/ ’ ^°  < VRI  < VRIO ^ 


( C-22 ) 


(C-23) 


Based  on  the  foregoing  derivations,  the  generalized  lamp- flasher  regulator 
input  characteristics  are  shown  in  Figure  C-8. 


VR|,  REGULATOR  INPUT  VOLTAGE,  vdc 


FIGURE  C-8.  LAMP-FLASHER  REGULATOR  INPUT  CHARACTERISTICS 


C-10 


FIGURE  C-9.  HOUSEKEEPING  REGULATOR  INPUT  CHARACTERISTICS 


C-ll 
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and  since  regulators  are  located  in  the  same  module 


TH  = rR 


Hence,  when  the  lamp  is  off: 


V ' IHI  {VRI'  V 

where:  1^  = User  load  current,  amperes 

which  condition  is  illustrated  in  Figure  C-10. 

When  the  lamp  is  on : 


I 


UL 


IHI  *VRI’  'V  + IRI  ^ VRI ’ TR^ 


which  condition  is  shown  in  Figure  C-ll. 


while  the  average  user  load  current  while  the  lamp  is  flashing  is 


TUL  IHI  ^VRI’  TR^  + XRI  ^ VRI  ’ TP.^ 


which  condition  is  given  in  Figure  C-12. 


Vg(  REGULATOR  INPUT  VOLTAGE,  .dc 
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Since  the  Lamp  Flasher  Unit  (user  load)  is  in  series  with  the  user  load 
cable : 


IPCD  JUL 


'-29 ) 


where:  Ip^,  = Power  Conditioning  and  Distribution  Group  current,  amperes 


VPCD  ^ IPCD*  TR*  " VRI  ^rPCD*  V + LPCD  RPCD  2-sL) 


where : 

VpcD  = Power  Conditioning  and  Distribution  Group  input  v , It age,  VDC 
RpCD  = User  load  cable  resistance,  ohms 

The  characteristics  of  the  power  conditioning  group,  based  on  Eqs.  (C-29) 
and  (C-30),  are  shown  graphically  in  Figure  C-13. 


VOLTAGE,  vdc 
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ENERGY  STORAGE  UNIT  ELECTRICAL  CHARACTERISTICS 
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The  simplest  form  of  the  Energy  Storage  Unit  in  relation  to  the  rest  of 
the  power  system  is  shown  in  Figure  D-l  below. 
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FIGURE  D-l.  ENERGY  STORAGE  UNIT,  NO  BATTERY  CHARGER 

The  relationship  between  the  battery  voltage  and  the  bus  voltage  (or  the 
battery  voltage  reflected  to  the  bus)  is  given  by: 


VBM  " VB  + TB  RL 


where : 


V = Battery  terminal  voltage,  VDC 
B 

V_w  = Battery  voltage  reflected  to  the  bus,  VDC 
BM 


I3  = Battery  current,  amperes 


+ = charge  mode 
- = discharge  mode 
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The  Energy  Storage  Unit  input  characteristics  under  these  conditions  is 
shown  in  Figure  D-2. 

An  Energy  Storage  Unit  with  a battery  charger  and  a battery  discharge 
diode  is  shown  in  Figure  D-3. 


During  periods  of  battery  charging,  the  operating  point  on  the  battery 
bus  is  shown  in  Figure  D-4. 

It  may  be  shown  that  the  operating  point  current  is  given  by: 


V - V 
CH00  BBM 


(D-2) 


or,  after  rearranging: 
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VB  ^BM 
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VOLTAGE,  vdc 
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FIGURE  D-2.  ENERGY  STORAGE  UNIT  CHARACTERISTICS,  NO  BATTERY  CHARGER 
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FIGURE  D-U . BATTERY  BUS  OPERATING  POINT  DURING  BATTERY  CHARGING 
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VCH00  = VBBM  + IBBM  ZCHR  fD_3) 

where: 

Vchoq  = Charger  output  voltage  at  zero  current,  amperes 
VBBM  = Operating  point  voltage,  VDC 
IgBM  = Operating  point  current,  amperes 
ZCHR  = Charger  output  impedance,  ohms 

The  regulation  characteristics  of  the  battery  charger  are  shown  in 
Figure  D-5. 

Examining  Figure  D-5,  it  can  be  seen  that: 

VESI  = <t>  { VCH00}  (D-1+) 


where:  VESI  = Energy  Storage  Unit  input  voltage,  VDC 
and  referring  to  Eq.  (D-3). 


FIGURE  D-5.  BATTERY  CHARGER  REGULATION  CHARACTERISTICS 

D-6 


r 
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In  the  saturated  region,  when  (Vr,HI0  < VESI  £ vcHIBA^: 


V = V + I Z 

CHOOS  BBM  BBM  CHRS 


(D-5) 


V = V + V 

ESI  CHOOS  OHIO 


(D-6) 


where : 


^CHR  ' = Charger  output  impedance  in  saturated  region,  ohms 

VCHOOS  = charSer  output  voltage  at  zero  current  in  saturated  region,  VDC 

V T.  = Input  voltage  at  charger  turn-on,  VDC 

V Vi  = Charge  input  voltage  at  changeover  from  saturated  to  active 
operation,  VDC 


While  in  the  active  region  of  operation:  (V_„T  •>  V ); 

iLO-L  LnlbA 


VCHOOA  = ^BBM  + TBBM  “CHRA 


VESI  = * *VCHOOA* 


(D-7) 


(D-8) 


where: 


VCHOoa  = c^arSer  output  voltage  at  zero  current  in  active  region,  VDC 
^CHRA  = Charge  output  impedance  in  the  active  region,  ohms 

Examining  Figure  D-3,  it  can  be  seen  that  during  battery  discharge 


V = V - V 
VESI  BM  DIODE 


where:  ^pjODE  = ®a'ttery  discharge  diode  voltage  drop,  VDC 


and  since  the  battery  charger  is  a series  regulator: 


IESI  = IB 


(D-9) 


( D-10 ) 


